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ABSTRACT 


E 1  ectromagnet i c  interference  from  power  lines  is  one  of 
the  main  sources  of  man-made  interference  to  communications 
in  the  high  frequency  <HF)  radio  band.  Two  types  of  radio 
interference  generated  by  power  lines  are  gap-type  noise 
caused  by  electric  discharges  across  line  hardware  and 
corona  noise  caused  by  the  partial  breakdown  of  the  air  due 
to  the  high  electric  fields  around  transmission  line 
conductors.  Using  original  data,  this  research  has  devel¬ 
oped  a  parameter  based  model  of  gap-type  and  corona  noise 
that  allows  the  fundamental  noise  mechanisms  to  be 
mathematically  or  physically  simulated.  An  expression  for 
the  power  spectral  density  (PSD)  of  gap-type  noise  and 
corona  is  derived.  The  energy  detection  problem  is 
formulated,  and  using  analytical  results  based  on  the  Hall 
model  for  radio  noise,  a  robust  energy  detection  receiver  is 
developed.  Tests  of  this  receiver  using  actual  and 
simulated  data  are  described. 
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I.  INTRODUCTION 


A.  PURPOSE 

Observations  of  radio  i nterference  at  high  ■frequency 
(HF)  receiver  sites  have  indicated  that  existing  models  of 
man-made  radio  noise  are  inadequate  to  describe  the  observed 
time-  and  frequency-domain  behavior  of  the  noise.  The  data 
also  shows  that  quite  often  a  specific  noise  source  pre¬ 
dominates  in  a  particular  location.  These  observations  led 
to  the  objectives  of  this  dissertation! 

1.  to  develop  specific  models  for  certain  man-made 
noise  sources  in  the  HF  radio  band  and 

2.  to  apply  this  model  to  the  analysis  of  energy 
detection  receivers.* 

B .  BACKGROUND 

Studies  conducted  over  the  past  few  years  have  indicated 
that  one  of  the  primary  sources  of  man-made  radio  noise  in 
the  HF  band  are  alternating  current  (AC)  transmission  and 
distribution  lines  CRefs.  1,23.  Two  of  the  primary  sources 
of  power  line  noise  are  gap  noise,  also  known  as  micro- 
sparking,  and  corona.  Gap  noise  is  caused  by  a  sparking 
process  between  hardware  points  on  utility  poles  and  corona 
is  caused  by  the  partial  breakdown  of  air  due  to  high 
electric  fields  around  high  voltage  conductors.  Both  of 
these  noise  types  are  non-Gaussian  noise  processes. 


The  majority  of  the  modal s  of  man-mada  and  atmospheric 
noise  to  data  have  assumed  the  high  amplitude  impulses 
driving  the  receiver  have  either  a  Poisson  arrival  rate  or  a 
variation  of  Poisson  arrivals.  This  assumption  has  allowed 
the  derivation  of  -first  order  envelope  statistics:  the 
amplitude  probability  distribution  (APD)  -function  and  the 
density  function  of  the  phase  of  received  noise  for  various 
amplitude  distributions  of  the  driving  impulses.  Observa¬ 
tions  of  man-made  noise  from  power  lines,  however,  has 
indicated  that  the  assumption  of  interpulse  independence 
(Poisson  arrival  times)  is  not  valid.  There  exists  a 
definite  time  domain  correlation  of  the  impulses  driving  the 
receiver.  This  is  due  to  the  underlying  deterministic 
mechanism  of  the  fundamental  frequency  of  the  power  line 
voltage.  Consequently,  this  research  has  concentrated  on 
the  statistics  of  the  impulse  arrivals  which  are  manifest  in 
the  autocorrelation  function  and  the  spectrum  of  the 
observed  noise  process.  Accordingly,  the  models  developed 
are  specific  rather  than  general  and  are  more  suited  to 
source  identification  by  spectrum  analysis  and  robustness 
evaluation  of  systems  rather  than  to  generic  specification 
of  optimum  receiver  structures. 

The  performance  of  energy  detection  receivers  for 
stochastic  signals  in  non-Gaussian  noise  has  received 
relatively  little  study.  Actual  implementation  of  algo¬ 
rithms  for  signal  detection  has  been  based  on  heuristics. 


The  advent  o F  digital  signal  processing  techniques  that 
allow  -For  complex  post-detection  algorithms  suggests  a 
care-fu  1 ,  stat i st i ca  1  analysis  of  the  noise  performance  of 
these  systems  may  help  develop  more  eFFicient  and  robust 
receivers. 


C.  CONTENT 

A  brieF  outline  oF  the  rest  oF  the  thesis  will  now  be 
presented.  Chapter  II  develops  a  general  model  For  bandpass 
impulsive  phenomena  using  a  Filtered  impulse  model  and 
complex  envelope  theory.  Nell  known  empirical  models  and 
physical  models  oF  radio  noise  are  presented  and  their 
relationship  to  the  generalized  model  discussed.  Particular 
attention  is  paid  to  the  Hall  model  For  atmospheric  radio 
noise  CReF.  33.  Some  useFul  extensions  to  it  are  derived 
that  will  be  used  For  simulations  in  Chapter  V. 

In  Chapter  III,  Field  observations  From  sources  oF  gap 


l _ 


noise  interFerence  are  presented  and  analyzed.  Three  cases 
oF  actual  interFerence  are  used  and  a  probabi 1 ist ic  model 
based  on  a  statistical  analysis  oF  the  data  and  the  Fi ltered 
impulse  model  developed  in  Chapter  II  is  speciFied.  The 
power  spectral  density  (PSD)  oF  gap  type  noise  as  predicted 
by  the  model  is  derived  and  compared  to  the  data. 

Using  an  actual  case  oF  corona  noise  From  a  500  kV  power 
line,  the  same  type  oF  analysis  is  carried  out  For  corona 
noise.  The  corona  noise  model  is  based  on  the  same  Fi ltered 


.  -  .v 
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impulse  -Framework  as  the  gap  noise  model  o-F  Chapter  III 
however,  the  spec  i-Fi  cat  ions  di-F-Fer  substantially. 

Using  the  Hall  model  -For  atmospheric  noise  discussed  in 
Chapter  II,  a  locally  optimum  receiver  -For  detection  of 
unknown  signals  in  HF  atmospheric  noise  is  derived.  A 
practical  modification  to  this  receiver  shows  it  to  be  an 
adaptive  limiter.  The  performance  of  this  receiver  is  then 
examined  in  simulated  man-made  noise  and  in  recorded  HF 
signal  and  noise  data. 

Appendix  A  presents  a  brief  description  of  the 
instrumentation  that  was  used  to  collect  much  of  the  data 
for  this  thesis.  The  dissertation  concludes  with  a  brief 
summary  of  results  and  some  suggestions  for  future  research. 


1 1 


jvaT.vT a." 
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INTRODUCTION 


In  -this  chapter  a  general  set  oF  spec iFi cat ions  -For  an 
impulsive  noise  interFerence  process  Mill  be  developed  and 
the  results  of  previous  work  in  the  Field  explained  in  terms 
oF  the  spec iFi cat ions.  This  generalized  scenario  will  be 
used  in  later  chapters  to  speciFy  complicated  sources  oF 
man-made  radio  interFerence .  In  attempting  to  describe 
atmospheric  radio  noise,  two  general  types  oF  models  have 
been  developed:  empirical  models  designed  to  Fit  First  order 
statistical  data,  and  physical  models  directly  related  to 
the  underlying  physical  mechanisms.  The  Hall  empirical 
model  For  atmospheric  radio  noise  has  been  shown  to  Fit 
atmospheric  noise  data  very  well  and  will  be  used  in  this 
dissertation.  Some  extensions  to  the  Hall  model  will  be 
developed  in  this  chapter. 

B.  GENERALIZED  SCENARIO 

In  order  to  provide  a  Framework  For  the  discussion  oF 
the  impulsive  noise  models,  a  general  interFerence  scenario 
For  impulsive  noise  will  be  described  CReF.  4,53.  A  typical 
interFerence  scenario  consists  oF  the  Following  elements: 

1.  a  source  oF  interFerence , 

2.  a  transmission  medium  to  the  receiver  and 
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3.  the  receiver  where  the  interFerence  man i Feats 
itselF. 

For  this  study  the  no is®  process  at  th®  r®c®iv®r  is 
modeled  as  th®  sum  oF  a  high  density  (in  tim®) ,  low  amp¬ 
litude  Gaussian  component  and  a  low  density,  high  amplitude 
impulsive  term. 

These  elements  are  shown  below,  where  e'(t)  in  Fig.  1  is 
the  impulsive  interference,  L'(w)  is  the  Frequency  response 
oF  the  transmission  medium,  z'(t>  is  white  Gaussian  noise 
and  H '  <u»)  is  the  combined  response  oF  the  RF  and  IF  Filters 
oF  the  receiver. 


z ' <t> 


e ' (t ) 


n ' (t) 


pw  M  U IJJ  *i.  IJWWfWWWWW  W L  ?  a  'A  W L'J "  V  M  JJ  *  U 


H  .  ’  .  l  H'J^IFilP.'PJ 


\  \ 


L^* 


The  source  of  impulsive  noise  is  specified  by  the 
equation 


N(t) 

e'(t>  -  E  e,  4(t-t,  ) 
i-1 


(2-1  ) 


where  e«  is  the  real  amplitude  of  the  ith  pulse  and  N(t)  is 
a  unit  counting  process  that  generates  arrival  times,  the 
ti 's.  Using  this  representation,  the  impulsive  interference 
source  is  described  by  the  probability  density  function 
(PDF)  of  each  et  ,  given  by  pt  (e>,  and  the  impulse  arrival 
times  generated  by  the  unit  counting  process,  N(t) . 

In  order  to  simplify  later  analysis  the  interference 
scenario  will  be  expressed  in  terms  of  its  complex  envelope. 
The  interference  process  at  the  receiver  detector  is  n'(t) 
and  will  be  replaced  by  its  complex  envelope  equivalent, 
n(t) ,  where 


jw*  t 

n'(t)  ■  ReCn(t)e  3 


and 


n«  (t)  +•  jn,  <t>  . 


(2-2) 


(2-3) 


n(t) 
t 

The  term  n(t)e  is  known  as  the  analytic  signal 
representat ion  of  n'(t),  and  n(t)  is  the  complex  envelope 
representation  of  n'(t).  The  reference  frequency  for  the 
complex  envelope  representat i on  is  wa  .  When  n'(t)  is  a 
bandpass  process,  ne (t)  and  n« (t)  are  the  lowpass  inphase 
and  quadrature  modulation  components  of  n(t)  respect¬ 
ively.!  Ref.  6:p.  753 


c 
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Specifying  a  filtered  impulse  model  in  complex  envelope 
form  requires  deriving  an  expression  for  the  driving  impulse 
function  given  in  Eqn.  2-1.  Two  equivalent  approaches  may 
be  taken  and  they  are  both  outlined  below.  In  the  first 
approach  a  bandpass  impulse  may  be  postulated  where  the 
spectrum  of  the  impulse  is  assumed  flat  from  u0  -  W  to  (u0  + 
W  where  ui0  is  the  reference  frequency  and  W  is  an  arbitrary 
bandwidth.  For  this  case  the  analytic  signal  representation 
is  easily  found  in  the  frequency  domain.  Using  the 
frequency  domain  definition,  the  analytic  signal  associated 
with  e'(t)  in  Eqn.  2-1  is 


e, (t)  =  (1/n) 


u> 

0  v 

E'  (title 


(2-4) 


Inserting  the  definition  of  the  bandpass  impulse  defined 
above  and  using  Eqn.  2-1 

N(t)  jot,  (t-t,  ) 

e,  (t)  =  E  (2/n )  e  sin(W(t-t, ) ) / (t-t,  > .  (2-5) 

i-1 

The  complex  envelope  of  e'(t)  is  defined  as 


-  ju«  t 


e(t)  =  e*  (tie 


(2-6) 


N (t )  -j0, 

e (t)  *  E  (2/n)  e  sin(W(t-t, ) ) / (t-t, )  (2-7) 

i  =  l 

where  otot,  =0,.  As  W  becomes  large,  e(t)  will  asympto¬ 
tically  approach  CRef.  71 


e(t)  ®  E  2e  6  (t-t,  ) 

i«l 


(2-8) 


The  ifeond  approach  [Ref.  83  is  to  consider  •the  impulse 
as  ideal  and  constant  in  -the  •frequency  domain.  The  analytic 
signal  representation  of  e'(t)  is  'then  directly  determined 
from  -the  •time  domain  definition  of  an  analytic  signal.  The 
analytic  signal  is  defined  as 


e«  (t>  «  e'(t)  je'C-t) 


(2-9) 


where  e'(t)  is  the  Hilbert  transform  of  e'(t).  The  analytic 
signal  representation  of  Eqn.  2-1  using  the  above  defin¬ 
itions  is 


N(t) 

e»  (t)  »  E  4(t-t,  )  +  j/(n(t-t,  )) 
i*l 


(2-10) 


since  1/nt  is  the  Hilbert  transform  of  the  delta  function. 
The  complex  envelope  of  e'(t)  is  then 


NCt)  -j0t  ~jw*t 

>(t)  3,1  E  &  (t-t»  )e  +  je  /(n(t-ti  >) 

1  =  1 


(2-11 ) 


using  Eqn.  2-6. 

Although  apparently  different,  when  the  complex  impulse 
trains  described  by  Eqns.  2-8  and  2-11  are  convolved  with  a 
complex  filter  impulse  response,  the  filter  output  is 
identical  for  both  representations.  If  h(t)  is  the  complex 
envelope  of  the  filter  impulse  response,  then  the  filter 


output  with  the  input  given  by  Eqn.  2-8  ii 


r 


E  © 

i  =  l 


h (t-t,  ) 


(.2-13) 


where  complex  convolution  is  defined  as 


h(t)  *  e(t)  =  (1/2) 


h  (  t  )  e  ( t  — r)dT. 


(2-13) 


The  -filter  output  with  the  input  defined  by  Eqn.  2-11  is 

N(t)  “j9i  -joiot 

(1/2)  E  h ( t-ti ) e  +  j©  h*  (t-t,  )  (2-14) 

i  =  l 

where  the  convolution  definition  of  the  Hilbert  transform 
was  used  CRef.  6:p.  693.  Using  Eqn.  2-6,  the  same  result  as 
derived  in  Eqn.  2-12  is  then  found.  Due  to  its  simplicity 
and  for  ease  of  use  in  programming  simulations,  Eqn.  2-8 
'will  be  used  to  represent  the  driving  impulse  function  for 
the  rest  of  this  dissertation. 

One  further  point  to  mention  concerns  the  probability 
distribution  of  8,  in  Eqn.  2-8.  For  the  processes  we  con¬ 
sider,  ti  will  have  random  arrival  times  and,  as  can  be  seen 
in  Eqn.  2-7,  8,  is  formed  by  multiplying  t,  times  the 
reference  frequency,  u0 .  Consequently,  8,  is  the  phase  of 
the  impulse  time  with  respect  to  the  reference  frequency. 
It  can  be  assumed  to  be  uniformly  distributed  over  0  to  2n 
when  the  t,  's  have  a  probab l 1 l st l c  interarrival  distribution 
and  Wo  is  much  greater  than  the  inverse  of  the  interarrival 
times  CRef.  9:  pp.  279,  101. 
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In  an  actual  interference  situation  the  received 
impulses  are  -filtered  by  the  radiating  antenna,  spreading 
losses,  atmospheric  attenuation,  receiving  antenna  and 
cabling  losses.  These  terms  will  be  lumped  together  as  an 
equivalent  filter,  L'(u).  Examples  of  this  filtering  term 
for  various  interference  sources  have  been  recorded  by  many 
researchers  CRef.  11,121.  The  filter  character  1 st 1 c  is 
specified  in  the  frequency  domain  as  L'(w)  since  that  is 
where  it  is  most  easily  observed. 

An  example  of  this  type  of  filtering  is  shown  in  Fig.  2 
where  the  frequency  spectrum  of  an  impulsive  noise  source 
due  to  a  gap  discharge  is  shown  from  0  to  200  MHz.  The 
straight  lines  in  the  3-Axis  view  are  due  to  stations  in  the 
high  frequency  (HF)  band  from  2  to  30  MHz  and  the  FM  band 
from  88  to  106  MHz.  The  continuous  envelope  seen  in  the 
upper  view  is  the  magnitude  of  the  frequency  response  of 
L'Cw).  The  solid  line  in  the  upper  picture  is  the  noise 
floor  of  the  spectrum  analyzer  at  -100  dBm. 

The  complex  envelope  of  the  impulse  response  of  the 
attenuation  term,  L'(u),  is  given  in  the  frequency  domain  by 
L(ti>)  =  L»CL'(-u)-(Do)  +  L'"(  ui+iJo  3  C  2  -  1  5  ) 

where  L*  denotes  the  lowpass  part  of  the  quantity  in 
brackets . 

A  complex  Gaussian  component  z(t)  will  also  exist  in  the 
observed  noise  process  and  is  due  to  the  combination  of: 


1.  the  thermal  noise  in  the  receiver,  and 


2.  the  combined  sum  of  many  low  level  atmospheric  or 
man-made  impulsive  sources. 

This  term  will  be  modeled  by  independent  zero  mean  Gaussian 
noise  processes  -for  the  inphase  and  quadrature  terms  with  a 
power  spectral  density  (PSD)  in  each  component  of  N0/2.  The 
distinction  between  the  Gaussian  noise  term  and  impulsive 
noise  term  is  that  many  receiver  responses  overlap  -for  the 
high  density,  low  amplitude  Gaussian  case,  satisfying  the 
condition  for  the  Central  Limit  Theorem.  For  the  impulsive 
noise  the  receiver  response  to  each  impulse  is  discernible 
such  that  the  probability  of  more  than  one  or  two  pulses 
overlapping  is  negligible.  Common  sources  of  impulsive 
noise  are  man-made  noise  due  to  power  lines  and  atmospheric 
noise  due  to  lightning. 

The  time  domain  input  to  the  receiver  is  now  modeled  by 
the  expression 


n  ( t ) 


(h(t)/2)  *  C  z (t ) 


N ( t )  -j8, 1  jw(t-t, 

+  E  2et  e  e 

i  =  1 


L  ( u) )  dbi  1 


(2-16) 


where 


h(t)  =  he(t)  jh,  (t). 

The  noise  power  bandwidth  of  H(u»)  will  be  designated  B*  r  t  . 
In  most  situations  L'(u>)  will  effectively  be  constant  in 
comparison  to  the  narrowband  receiver  filter  H(u).  This 
makes  L(uj)  a  constant  that  depends  only  on  w0  . 
Additionally,  g(t)  will  be  defined  by 
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g  (t ) 


(h(t/2))  *  z(t) 


(2-17) 


and  is  -the  low  level  Gaussian  noise  at  the  detector.  The 
power  in  the  inphase  and  quadrature  components  of  the 
Gaussian  noise  will  be  NoB«**/2.  The  complex  lowpass  noise 
process  can  now  be  written  as 

N(t)  -j9, 

n(t)  **  g(t>  E  a,  e  h(t-t,  )  (2-18) 

i»l 


where 

at  =»  e,L '((!)«).  (2-19) 
The  parameter  a,  is  then  a  scenario  dependent  parameter  that 
is  a  -function  o-f  the  reference  frequency  and  the  impulse 
amplitude.  With  the  above  simplifications,  Fig.  1  can  now 
be  expressed  in  complex  envelope  form  as  shown  in  Fig.  3. 


z  (t) 


Figure  3.  Complex  Envelope  of  Interference  Scenario 

The  terms  from  Fig.  1  have  been  replaced  by  their  equivalent 
complex  envelope  representations  and  the  prime  is  dropped  to 


Some  -further  simplifications  will  now  be  considered 
based  on  the  characteristics  of  the  H'(cu),  the  bandpass 
filter.  If  the  impulse  response  hCt)  of  the  f 1  Iter  is  real, 
H'(uj)  has  conjugate  symmetry  about  the  origin  in  the  sense 
that  H'*(ui)  =  H'(-tj).  Furthermore,  if  H'(w)  is  sym-metrical 
in  the  same  way  about  the  reference  frequency,  ,  then 

h, (t)  will  equal  0. 

The  envelope  squared  of  the  process  defined  in  Eqn.  2-18 
1  s 

2  2  2 

E  (t)  =  nc  (t)  +  n,  (t).  (2-20) 

If  the  Gaussian  noise  term  is  assumed  to  equal  0  then 
2  NCT)  NCT) 

E  (t)  =  E  E  a,  a,  hc  (t-t,  )  hc  Ct-tj  )cos(8,  -0,  )  + 

1=1  j-1 

NCT)  NCT) 

E  E  a,  aj  h,  (t-t,  )ht  (t-t,  )  cos  (9,-0j).  (2-21) 

i=l  j-1 

If  NCT)  defines  a  low  density  counting  process  such  that 
there  is  negligible  overlap  between  subsequent  pulses  and 
h(t)  is  the  complex  envelope  of  the  impulse  response  of  a 
bandpass  symmetrical  filter,  then  the  envelope  of  n(t)  can 
be  simply  written  as 

N(t) 

E (t)  =  E  a, h ( t-t, )  (2-22) 

i  =  l 

where  h(t)  is  equal  to  hc  (t) . 

At  this  point  of  the  analysis  two  avenues  may  be 
pursued.  If  there  is  a  time  dependency  between  pulses  (non- 


exponent i a  1 1 y  distributed  interarrival  times),  then  tne 
correlation  function  and  spectrum  of  the  process  can  be 
examined.  If  the  amplitude  probability  distribution  (APD) 
o f  the  envelope  is  the  desired  result,  then  the  char¬ 
acteristic  function  of  the  envelope  is  most  useful  .  The  APD 
function  is  commonly  used  in  radio  noise  research  and  is 
equal  to  one  minus  the  cumulative  distribution  function 
(CDF).  In  general  these  two  approaches  are  mutually 
exclusive  since  the  assumption  of  time  dependency  between 
the  pulses  makes  the  envelope  APD  function  difficult  to 
calculate.  If  the  impulse  arrival  times  form  a  Poisson 
point  process  then  the  envelope  APD  may  be  calculated: 
however,  the  spectrum  is  constant.  The  remainder  of  this 
chapter  will  follow  the  envelope  APD  approach.  Results 
obtained  by  previous  researchers  will  be  developed  using  the 
general  noise  model.  Chapters  III  and  IV  will  explain  the 
other  avenue  and  look  at  the  spectrum  of  complicated  man¬ 
made  noise  processes. 

C.  EMPIRICAL  MODELS 

Empirical  models  of  atmospheric  noise  have  been 
developed  to  provide  a  mathematical  expression  for  the 
first-order  statistics  of  the  envelope  of  the  received 
waveform.  In  particular,  the  APD  has  been  emphasized.  This 
type  of  model  attempts  to  construct  a  mathematical  express¬ 
ion  that  fits  observed  data  without  regard  for  the  physics 


of  "the  i nterference  scenario 


Empirical  models  have  "the 


advantage  that  the  resulting  expressions  are  much  simpler 
than  those  obtained  -Prom  the  filtered  impulse  models.  One 
disadvantage  of  this  class  of  model  is  that  since  only  the 
first  order  statistics  of  the  noise  are  considered  in 
developing  and  fitting  the  model,  the  higher  order  sta¬ 
tistics  may  not  match  the  data  well. 
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p(r)  =  2ore 

r  2  0 

(2-23) 

=  0 

r  <  0 

and  the  APD 

or 

exceedance  probability 

2 

-ar, 

Pr  (r 

>  r0  )  =  1  -  P(r0  )  =  e 

r*  2  0 

(2-24) 

=  1 

r„  <  0 

How  well  the  single  parameter  Rayleigh  distribution  fits 
observed  atmospheric  noise  data  can  be  seen  in  Fig.  4.  The 
data  points  are  exceedance  probabi 1 i t ies  plotted  in  decibels 
(dB)  above  the  root-mean-squared  (RMS)  value  of  the  received 
envelope.  The  data  was  measured  on  an  ARN-2  receiver  at  a 
center  frequency  of  10  MHz  at  Boulder,  Colorado  [Ref.  133. 
The  high  probability,  low  amplitude  portions  of  the  observed 
data  approaches  the  slope  of  the  Rayleigh  distribution  curve 
for  a  equal  to  12.  (a  is  defined  in  Eqn.  2-24)  This  value 
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Figure  4.  Atmospheric  No is*  and  Rayleigh  Distribution 


oF  a  was  chosen  to  Fit  the  high  probability  portion  oF  the 
data.  That  the  data  approaches  a  Rayleigh  distribution  is 
due  to  the  e-FFect  oF  the  many  small  overlapping  impulses 
occurring  at  low  amplitudes  whose  quadrature  Gaussian 
components  will  have  a  Rayleigh  envelope.  However,  at  the 
low  probability,  high  amplitude  part  oF  the  curve  the 
Rayleigh  distribution  predicts  Far  Fewer  values  than 
actually  occur.  This  is  the  distinguishing  Factor  oF  the 
envelope  distribution  oF  impulsive  noise:  that  its 
distribution  contains  a  higher  probability  oF  high  amplitude 
terms  than  is  predicted  by  the  Rayleigh  distribution.  The 
second  curve  in  Fig.  4  where  a  is  equal  to  .5  is  the 
Rayleigh  distribution  with  the  same  power  as  the  observed 
data.  It  can  be  seen  that  this  curve  under  estimates  the 
high  amplitude  and  over  estimates  the  low  amplitude  values 
oF  the  observed  distribution.  It  should  be  noted  that  this 
type  oF  plot  emphasizes  the  low  probability  portion  oF 
distribution  since  that  is  where  the  deviation  From  the 
expected  behavior  occurs. 

In  an  attempt  to  correct  the  poor  Fit  oF  the  Rayleigh 
distribution  at  high  amplitude  levels,  Likhter  CReF.  143 
proposed  a  combination  oF  two  Rayleigh  distributions: 

2  2 
-or»  -8r» 

Pr(r  >r0)=<l-c)e  +  ce  .  (2-25) 
This  Formula  has  been  shown  to  agree  poorly  with  practical 
results  CReF.  153. 
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Spaulding,  Roubique  and  Crichlow  CReF.  161  combined  the 
Rayleigh  distribution  with  a  "power"  Rayleigh  distribution 
to  obtain  a  distribution  that  Fit  very  well  For  atmospheric 
noise  over  a  wide  range  of  receiver  bandwidths.  This  APD  is 

2 

-otr. 

Pr  <  r  >  r«  )  =  e  r®  <  fl 

2  1/s 
-(or®  ) 

Pr  (r  >  r®  )  =  e  r®  2  fl  (2-26) 

where  o,  fl  and  s  are  determined  From  measured  statistical 
parameters  oF  the  noise. 

Horner  and  Harwood  CReF.  173  proposed  the  two  parameter 
log-normal  dii.  ribution  and  Found  it  gave  a  satisFactory  Fit 
to  radio  noise  data  in  the  VLF  band.  The  log-normal  was 
chosen  since  it  has  a  more  impulsive  tail  and  Fits  the  high 
amplitude,  low  probability  data  better  than  the  Rayleigh 
distribution.  The  PDF  For  the  log-normal  distribution  is 

2  2 

1  -(ln(r)-o)  /2o 

p<r)  »  -  e  (2-27) 

1/2 

r<7  (2n  ) 

Two  examples  oF  the  APD  oF  the  log-normal  distribution  are 
shown  in  Fig.  5  where  the  curve  with  a  equal  to  -.24  and  a* 
equal  to  .24  is  the  log-normal  distribution  with  the  same 
mean  and  mean  square  parameters  as  the  curve  with  a  equal  to 
.5  in  Fig.  4.  Comparing  the  two  curves  the  log-normal 
distribution  provides  an  intrinsically  better  Fit  to 
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impulsive  noise  date  at  both  the  high  and  low  amplitude 
values.  When  the  log-normal  density  has  both  parameters 
chosen  to  match  the  mean  and  mean  square  values  of  the 
observed  data  an  excellent  fit  to  this  particular  data  is 
obtained.  This  is  seen  in  the  curve  with  a  equal  to  -1.4 
and  <r*  equal  to  1.4. 

Another  model  that  has  been  successfully  used  to  fit  the 
observed  APD  of  atmospheric  radio  noise  is  Hall's  general¬ 
ized  "t"  distribution  CRef.  33.  This  model  is  unique  in 
that  it  is  not  of  the  filtered  impulse  type;  however,  in 
contrast  to  the  strictly  empirical  models  it  does  postulate 
a  random  process. 

One  of  the  problems  of  Gaussian  models  of  atmospheric 
noise  is  that  the  dynamic  range  of  the  model  is  less  than 
the  observed  dynamic  range  of  the  noise.  To  achieve  a 
greater  dynamic  range  starting  with  a  Gaussian  process,  Hall 
proposed  a  model  which  considered  the  received  noise  to  be 
of  the  form 

n<t )  »  a<t)s(t)  (2-28) 
where  a<t)  is  a  slowly  varying  stationary  random  process  and 
s(t)  is  an  independent,  narrowband  Gaussian  process.  An 
analytically  tractable  distribution  was  selected  for  a(t) 
which  was  chosen  to  give  good  agreement  between  the  model 
and  measured  atmospheric  noise  data.  The  distribution  for 


a(t)  is: 


<m/2>  1  -(m/2a  a  ) 

p  ( a  )  ■  ■—  -  —  • 

m  m+1 

a  r (m/2)  | a |  (2-29) 

where  m  and  a  are  the  two  parameters  used  to  fit  the  model 
to  the  data.  The  distribution  of  s(t)  is  Gaussian; 


2  2 
<s  /2<r,  ) 


pCs) 


2  1/2 
<2nffi  > 


(2-30) 


Hall  calculated  the  density  function  of  n(t)  to  be 

0-1 


r<0/2) 


1 


p  <n) 


1/2  2  2  0/2 
T(  <0-1  )/2)  n  (n  Y  > 


(2-31 ) 


where 


1/2 

y  »  m-  (<r«  /<r )  , 


0  *  m  ■  1 , 

and  r  is  the  gamma  function.  For  the  case  where  ffi  =  <r , 
Eqn.  2-31  is  the  density  function  for  Student's  "t" 
distribution  which  is  the  basis  for  describing  the  density 
as  a  generalized  "t"  distribution. 

Hall  then  calculated  the  envelope  distribution  for  n(t) 
based  on  the  above  assumptions  and  obtained 


0-1 


p(r ) 


<0  -  1)  y 


(r 


2  (<0+l)/2> 
Y  > 


(2-32) 


with  the  phase  uniformly  distributed  between  0  and  2n . 


For  atmospheric  radio  noise,  values  of  8  in  the  range  2 


to  5  have  been  Found  to  give  excellent  agreement  with  the 
data  CReF.  18,193.  (Models  with  integer  values  oF  the 
parameter  8  will  be  abbreviated  to  Hal  12,  Hal  13,  etc.) 

Fig.  6  shows  the  Hal  13  Fit  to  the  same  atmospheric  noise 
data  used  in  Figs.  3  and  4.  An  excellent  Fit  to  this  data 
is  obtained.  The  Hall  model  has  some  disadvantages.  The 
higher  moments  only  exist  For  orders  greater  than  8-1  and 
the  parameters  8  and  y  must  be  determined  For  each 
interFerence  scenario  and  are  not  easily  related  to  the 
physical  source  oF  the  interFerence . 

D.  FILTERED  IMPULSE  NOISE  MODELS 

The  Filtered  impulse  models  diFFer  Fundamentally,  in 
concept  From  the  empirical  models  described  previously. 
This  class  oF  model  will  be  explained  in  terms  oF  the 
generalized  interFerence  scenario  presented  at  the  beginning 
oF  the  chapter.  The  Filtered  impulse  models  have  the 
advantage  oF  being  based  on  the  underlying  physics  oF  the 
process  but  suFFer  From  the  disadvantage  oF  being 
analytically  complex.  Interestingly,  some  oF  the  results 
obtained  From  the  Filtered  impulse  models  have  retroactively 
justiFied  the  expression'.  that  were  derived  From  the 
empirical  models. 

To  determine  the  APD  oF  the  envelope  oF  a  Filtered 
impulse  process,  the  joint  character istic  Function  oF  the 
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inphase  and  quadrature  components  is  determined. 


Since  the 


noise  is  narrowband  it  will  be  circularly  symmetrical  and 
the  2-D  Fourier  transform  will  become  a  Hankel  transform 
CRef.  20,213.  This  character i st i c  function  is  transformed 
to  polar  coordinates,  then  inverted  to  determine  the  APD . 
Starting  with  the  complex  noise  process  described  in  Eqn.  2- 
18,  the  joint  character i st i c  function  of  the  inphase  and 
quadrature  components  is 

j  (a),  nc  (t)  +  in,  n,  (t)  ) 

0<(i)i  |iiii  )  =  ECe  3.  (2-33) 

It  is  shown  in  Appendix  B  that  the  character l st l c  function 
of  the  envelope  due  to  the  Gaussian  term  is 

2 

-No  B  ( to,  )/2 

0a  (u)r  )  =  e  (2-34) 

and  the  character i st i c  function  of  the  envelope  due  to  the 
impulsive  term  is 

00  T 

X  p  ( a )  C  Jo  ( w,  ahc  (t  )  )  -  lldTda  (2-35) 

0  0 

0i  ( ui,  )  =  e 

where  p(a)  is  the  impulse  amplitude  density  function,  X,  the 

Poisson  rate  function,  he  (t),  the  inphase  component  of  the 

filter  impulse  response  and  J0  is  the  ordinary  Bessel 
function  of  the  first  kind.  The  density  function  of  the 

envelope  is  then  the  inverse  Hankel  transform  given  by 

00 

p(r)  =  ru),  Jo(rii>r)0*((i»r)0i(u)r)du>r.  (2-36) 

v 
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2.  Uni -form  spatial  distribution  oF  impulse  sources 


3.  Inverse  distance  law  -For  received  Field  strength 

4.  Arbitrary  receiver  envelope  response 

and  showed  that  the  APD  oF  the  envelope  Function  was  approx- 
i mated  by 

K 

Pr(r  >  r.  )  »  -  (2-37) 


2  2  1/2 


(r.  +  K  ) 


where 


T 


K 


h (t )  dt 


0 


T  is  the  observation  interval,  h(t>  the  Fi Iter  response,  and 
c  is  a  scenario  dependent  constant.-  This  result  is  sig- 
niFicant  because  it  has  the  same  Form  as  the  APD  oF  the 
Hal  12  model  and  physically  justiFies  what  had  heretoFore 
been  an  empirical  model. 

Middleton  [ReF.  24]  included  the  Gaussian  noise  term  in 
Eqn.  2-36  and  expanded  the  characteristic  Function  without 
taking  the  expectation.  Then  by  inverting  the  char¬ 
acteristic  Function  term  by  term,  a  canonical  Form  oF  the 
envelope  density  Function  can  be  obtained.  The  exceedance 
probability  was  then  shown  to  be 
Pr(r  >  r»  )  = 

m 

2  2  •  ( - 1 )  m  mo  mo 

-  r»  Cl  -  r»  E  -  A.  T  (  1  +  - ) ,  F,  (  1  -  - ;  2 ;  r<>  )  ] 

m=0  m !  2  2 


(2-38) 


where  i  F,  is  a  conFluent  hypergeometr ic  function,  and  a  and 
A.  art  'two  parameters  that  are  determined  by  the  source  dis¬ 
tributions  and  the  propagation  law. 


E.  HALL  MODEL  EXTENSIONS 


For  later  work  in  this  dissertation  a  model  For 
atmospheric  noise  statistics  is  required.  This  section  will 
Further  describe  the  Hall  model  as  applied  to  HF  atmospheric 
noise  and  will  derive  some  extensions  to  it  that  will  be 
used  to  digitally  simulate  atmospheric  radio  noise. 

The  Hall  model  For  the  envelope  statistics  oF  radio 


noise  is 


0-1 


r 


p<r>  =  <0  -  1)  y  - 

2  2  (C0-l)/2) 

<r  +  y  >  (2-39) 

and  was  chosen  For  a  number  oF  reasons.  The  most  important 

arei  the  Hall  model  is  analytically  tractable  and  shows  a 

good  Fit  to  the  data.  It  can  be  obtained  For  certain  values 

oF  the  parameter  0,  starting  From  the  physically  based 

Filtered  impulse  models.  This  was  shown  by  both  Giordano 

CReF.  233  and  Middleton  CReF.  243.  Additionally,  SchonhoFF 

CReF.  181  has  related  the  Hall  models  to  CCIR  Report  332 

CReF.  253,  the  standard  reFerence  oF  First  order  statistics 

oF  atmospheric  noise.  In  the  CCIR  report  the  APD  curves  are 

shown  as  a  Function  oF  V0 .  The  parameter,  V*  ,  is  a  commonly 

used  measure  oF  the  impulsiveness  oF  atmospheric  radio  noise 

and  is  deFined  as 


*  -  ■  -  ■ 


■V-V-V. 
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2  1/2 
ECr  1 

VD  =  201og  - 

ECr]  (2-40) 

where  r  is  the  amplitude  oF  the  received  envelope.  For 

quadrature  Gaussian  noise  VD  =  1.05  dB  and  as  the  noise 

envelope  becomes  more  impulsive,  V„  will  increase. 

One  problem  in  relating  the  Hall  model  to  VD  is  that  the 
moments  o-F  the  Hall  model  only  exist  For  orders  greater  than 
8-1.  This  means  that  For  the  Hall2  and  Hall3  models,  Vo  is 
undeFined  unless  the  model  is  modiFied.  Two  methods  have 
been  proposed  to  do  this.  Hall's  method  CReF.  31  adds 
another  term  to  the  envelope  density  Function  and  multiplies 
by  a  negative  exponential  term  to  reduce  the  tail  oF  the 
density  Function  and  give  Finite  moments.  The  advantage  oF 
this  method  is  that  the  density  Function  exists  From  0  to 
inFinity;  however,  the  resulting  density  Function  is 
complicated.  The  method  proposed  by  SchonhoFF  CReF.  13] 
truncates  the  Hal  12  and  Hal  13  distributions  above  a  set 
level,  TP  ,  and  renormalizes  the  density  Function  by  picking 
k  such  that 

Tp 

k  p(r)dr  =  1.  C  2  —  <4  1  'i 

0 

The  level,  TP  ,  is  chosen  to  set  a  desired  Vo  For  the 
distribution.  The  constant  k  will  be  a  Function  oF  D  where 

2  1/2 

=  (1  +  (Tp /y)  ) . 


0 


(2-42) 


The  existence  oF  the  -First  and  second  moments  -For  these  two 


distributions  is  then  assured.  Using  this  procedure 
SchonhoFF  generated  a  -Family  o-F  distributions  parameter  1  zed 
by  0  that  can  be  used  to  represent  a  wide  range  o-F  measured 
data.  By  matching  the  ratio  o-F  the  data  to  the  dis¬ 
tribution,  a  close  -Fit  to  the  APD  is  -Found. 

An  extension  to  the  basic  Hall  model  to  be  used  in  a 
later  chapter  is  the  density  -Function  o-F  a  constant 
amplitude  sinusoid  plus  the  assumed  noise.  This  density  is 
analogous  to  the  Rician  density  in  Gaussian  noise  theory. 
It  will  be  used  in  simulating  the  performance  o-F  energy 
detection  receivers  -For  the  signal  plus  noise  case.  It  was 
derived  in  CRe-F.  263  -For  the  Hal  12  model  and  it  is  derived 
-For  the  Hal  13  and  Hal  15  models  in  Appendix  C.  For  the  Hal  13 
model  the  density  -Function  is 


2  2  2  2 

<  r  +  A  +  y  )  2y  rD 

p  ( r )  =  -  (2-43) 

4  22  22  22  242 

<r  -  2A  r  +  2A  y  2+  2y  r  2+  A  +  y  ) (D  -  1) 


where  A  is  the  signal  amplitude,  y  is  a  Hal  13  noise 

parameter  and  D  is  de-Fined  in  Eqn.  2-42.  Likewise  the 

density  Function  For  the  Hal  15  model  is 

4  2  2 

2y  r  <  b  +  2a  ) 

p(r)  =  -  (2-44) 

2  2  5/2 

(a  -  b  ) 

where 


Fig.  7  is  an  example  of  the  aPD  for  a  sinusoid  plus  Hal  13 
noise  with  the  RMS  value  of  the  noise  set  to  1  and  the 
signal  to  noise  ratio  at  -26,  -6  and  14  dB.  For  the  -26  dB 
SNR  case  the  APD  is  very  similar  to  the  Hal  13  noise  only 
density  shown  in  Fig.  6. 

Commonly  used  functions  relating  to  the  Hall  models  are 
shown  in  Table  I.  These  are  the  envelope  density  function, 
the  1st  and  2nd  moments  of  the  envelope,  the  inverse  of  the 
CDF,  the  quadrature  marginal  density  function  and  the  Hall 
noise  plus  random  sinusoid  density.  The  parameter  D  in 
Table  I  is  defined  by  Eqn.  2-42  and  the  function  £C • D  is  the 
complete  elliptic  irrtagral  of  the  second  kind.  The  quad¬ 
rature  components  can  be  derived  by  transforming  the 
envelope  density  to  rectangular  co-ordinates  and  integrating 
with  respect  to  one  component.  The  inverse  of  the  envelope 
CDF  is  also  included  since  it  provides  an  easy  method  to 
generate  random  deviates. 

Also  useful  for  simulation  purposes  are  the  values  of 
the  Hall  parameters,  0  and  y ,  for  different  values  of  V0 
normalized  to  a  unit  root  mean  square  value.  The  truncation 
point  Tp  was  found  for  a  desired  V»  using  an  iterative  tech¬ 
nique.  These  parameters  are  shown  in  Table  II. 
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TABLE 


V*  Hall  No.  ECn  <t>3  ECn(t)l  T, 


2. 1 

5 

1 

.785 

- 

1  .01 

- 

3 

4 

1 

.708 

- 

.708 

- 

4 

3 

1 

.631 

11.66 

.420 

27.76 

N 

•\V. 

5 

3 

1 

.562 

25.90 

.364 

71  .05 

6 

3 

1 

.501 

67.71 

.321 

211.1 

8 

2 

1 

.398 

11.77 

.085 

137.6 

10 

2 

1 

.316 

17.08 

.058 

291.0 

\V.V 

\«V 

12 

2 

1 

.251 

24.25 

.041 

587. 1 

\V.‘- 

14 

2 

1 

.  199 

33.81 

.029 

1141.9 

S  '<4 

TABLE  II 


F .  SUMMARY 

Both  empirical  and  analytical  Filtered  impulse  models 
have  been  successful  in  explaining  the  First  order  envelope 
statistics  oF  bandpass  radio  noise  processes.  However,  in 
obtaining  these  results  the  assumption  has  to  be  made  that 
the  impulse  arrival  times  due  to  the  noise  process  Form  a 
Poisson  point  process.  In  general  this  will  not  be  true  For 
man-made  impulsive  noise  and  For  some  types  oF  atmospheric 
noise  CReF.  271.  In  most  cases  oF  man-made  noise  the 
underlying  deterministic  mechanism  oF  the  device  causing  the 
noise  will  modulate  or  regulate  the  impulse  distribution. 
The  Hall  model  For  atmospheric  noise  has  been  examined  in 
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III.  GAP  NOISE  MODEL 


A.  INTRODUCTION 

In  'this  chapter  a  specific  source  of  man-made  radio 
noise  will  be  modeled  and  specified  in  terms  oF  the  counting 
process,  N<t) ,  and  the  probability  density  function  (PDF)  oF 
the  impulses,  p<a) .  These  spec iFi cat ions  were  introduced  in 
the  previous  chapter  to  deFine  the  impulsive  noise  model. 
The  noise  source  that  will  be  examined  is  gap  type 
discharges  and  is  commonly  Found  on  electric  power  dis¬ 
tribution  and  transmission  lines.  Gap  noise  is  one  oF  the 
major  types  oF  i nterFerence  From  power  lines  and  is 
Frequently  observed  as  the  primary  i nterFerence  to  com¬ 
munication  systems  operating  in  the  high  Frequency  (HF> 
radio  band. 

B.  GAP  NOISE  MECHANISM 

At  least  two  mechanisms  have  been  Found  by  which  a  gap 
discharge  process  can  occur  on  a  power  line.  The  resistance 
in  the  line  insulators  can  be  degraded  allowing  current  to 
Flow  through  the  insulator  base  and  creating  a  potential 
gradient  across  any  gaps  or  deFects  in  the  insulator 
mounting  hardware.  A  second  way  in  which  a  potential  can  be 
created  across  an  air  gap  is  by  an  electro-stat ic  coupling 
oF  the  line  potential  to  isolated  hardware  on  the  pole.  In 


both  cases  the  potential  across  the  gap  is  discharged  by  the 


voltage  breakdown  oF  the  gap  and  the  resulting  rapid  current 
■Flow  or  spark.  This  process  generates  a  radio  Frequency 
(RF)  noise  impulse  with  spectral  components  extending  into 
the  hundreds  oF  MHz.  During  a  single  discharge,  the 
potential  across  the  gap  is  temporarily  diminished.  How¬ 
ever,  while  the  Fundamental  60-Hz  waveForm  is  still  greater 
than  a  threshold  voltage,  the  process  can  occur  again.  The 
spark  will  discharge  across  the  gap  repeatedly  until  the 
alternating  current  waveForm  drops  below  the  breakdown 
threshold  potential.  CReF.  28: p.  78,291 

The  spark  discharge  and  recharging  oF  the  gap  potential 
indicates  that  this  type  oF  process  is  regenerative  and  can 
be  modeled  as  a  renewal  process  where  the  renewal  points  are 
associated  with  the  sparks.  One  Feature  oF  gap  noise  that 
complicates  the  modeling  is  the  120-Hz  on-oFF-on  modulation 
imposed  on  the  renewal  points  by  the  alternating  current 
(AC)  waveForm.  One  way  to  account  For  this  eFFect  is  to 
consider  the  turn  on  times  as  another  renewal  process  driv¬ 
ing  the  spark  discharge  process.  This  type  oF  model  is 
known  as  a  branching  renewal  process  CReF.  301. 

The  noise  processes  that  were  used  to  develop  this  model 
were  short  term  stationary  in  that  the  statistics  oF  the 
process  did  not  change  over  the  observation  interval.  Not 
all  sources  oF  power  noise  are  stationary  even  over  a  short 
observation  interval.  Physical  eFFects  such  as  wind,  solar 
heating  and  varying  line  loads  can  act  to  make  certain  types 


o-f  power  line  noise  highly  variable  'from  observation  to 


observation.  Although  these  types  of  noise  Mere  observed 
during  data  collection,  the  data  used  to  develop  the  model 
was  taken  from  noise  sources  that  were  stationary  for  the 
length  of  the  ten  minute  data  records. 


C.  OBSERVED  TIME -DOHA IN  CHARACTERISTICS 

The  gap  noise  sources  for  this  study  were  observed  on 
utility  distribution  lines  in  the  vicinity  of  the  Naval 
Postgraduate  School  and  were  chosen  to  illustrate  parameters 
of  the  noise  model.  Fig.  8  is  a  typical  time-domain  observ¬ 
ation  of  a  gap  discharge  process  observed  at  3  MHz  and 
envelope  demodulated  with  a  10-kHz  Gaussian  bandpass  filter. 
The  important  characteristics  to  note  are: 

1.  the  process  has  an  on-off-on  modulation  at  a  120- 
Hz  rate  related  to  60-Hz  waveform  of  the  power 
line,  and 

2.  the  pulse  groups  that  result  from  the  modulation 
have  a  probabilistic  number  of  impulses  occurring 
in  each  group  and  a  probabilistic  interarrival 
time  between  pulses  in  a  group. 

In  the  first  pulse  group,  9  pulses  occur  with  varying 

amplitudes  and  interpulse  arrival  times.  In  the  second 

pulse  group  starting  approximately  8.33  msec  later,  8  pulses 

occur,  again  with  varying  amplitudes  and  interarrival  times. 

One  group  of  pulses  is  associated  with  either  the  positive 

or  negative  polarity  of  the  line  voltage  waveform  and  the 

other  group  with  the  opposite  polarity.  Identification  of 
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vs 


'the  polarity  is  impossible  without  physically  locating  the 


source . 

The  average  amplitude  of  the  impulses  will  be  one  of  the 
parameters  -for  the  model.  An  estimate  o-f  this  parameter  can 
be  made  -from  the  above  presentation .  It  is  important  to 
note  that  -for  this  highly  impulsive  type  o-f  noise  that  the 
observed  amplitude  is  a  -function  o-f  the  shape  and  width  of 
the  effective  bandpass  filter  and  the  detector  char¬ 
acter  ist  ics  . 

At  this  point  it  should  be  noted  that  the  model  can  be 
approached  at  two  levels  of  complexity.  In  the  simple 
version  the  data  from  both  the  positive  and  negative  phases 
is  considered  as  a  whole  and  averaged  over  the  two  phases  to 
determine  a  set  of  average  parameters.  In  the  more  complex, 
but  more  accurate  model,  two  sets  of  model  parameters  are 
determined  -  one  for  each  phase.  For  the  remainder  of  this 
work  the  more  complex  case  will  be  considered.  The  results 
for  the  simple  case  are  presented  in  Moose  and  O'Dwyer  [Ref. 
313  and  can  easily  be  determined  from  the  complex  case  by 
setting  the  parameters  equal  for  both  phases  and  simplifying 
the  resulting  equations. 

In  order  to  better  characterize  the  interpulse  arrival 
times,  which  will  be  used  to  determine  two  model  parameters, 
the  rising  raster  capability  described  in  the  instru¬ 
mentation  section  was  used  to  generate  the  display  of  a 
different  gap  noise  source  which  is  shown  in  Fig.  9.  In 
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■this  picture  the  amplitude  oF  each  individual  record  was 


reduced  and  thresholded  so  that  all  that  essentially  remains 
is  the  time  oF  impulse  arrival  data  For  28  time  records. 
The  time  base  oF  the  display  was  intentionally  synchronized 
to  the  power  line  to  Facilitate  taking  data  and  this 
accounts  For  the  regularity  oF  the  pulse  groups  From 
observation  to  observation.  In  this  view  the  interarrival 
times  For  56  pulse  groups  (28  groups  oF  one  polarity  and  28 
oF  the  opposite  polarity)  can  be  determined  along  with  the 
number  oF  pulses  in  each  oF  the  56  successive  pulse  groups. 
The  average  number  oF  pulses  per  group  considered  separately 
For  each  phase  will  also  be  used  as  estimates  o F  two  model 
parameters . 

Fig.  10  is  a  histogram  oF  the  distribution  of  the 
interarrival  times  between  the  observed  impulses  For  the 
same  noise  process  shown  in  Fig.  9.  To  use  this  data  to 
generate  parameters  For  the  noise  model,  the  histogram  will 
be  approximated  by  a  continuous  density  Function.  The  gamma 
density  Function  CReF.  32:p.  183 

r-1  -CXt) 

p(t)  =  X(Xt)  e  /(r-1)!  (3-1) 
was  chosen  because  it  showed  a  close  Fit  to  the  data ,  and 
its  character istic  Function,  which  will  be  used  in  later 
derivations,  was  particularly  simple.  The  characterist ic 
Function  oF  the  gamma  density  is 

-r 

0(i<u)  =  (1  -  -.uX)  .  (3-2) 
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The  r  and  X  parameters  are  simply  interpreted  in  terms 


oP  an  underlying  Poisson  process  as  the  time  to  the  rth 
point  oP  a  Poisson  process  of  intensity  X.  IP  p  is  the 
random  variable  assigned  to  the  observed  interarrival  times 


then  the  unbiased  estimate  of  the  mean  of  p  is 


est(p)  =  (1/N)  E  p, 
i  =  l 

and  the  unbiased  estimate  of  the  variance  is 

2  N  2  _  2 

est(aH  )  *  {  E  pi  -  N  estCp)  >/(N-l). 
1  =  1 


(3-3) 


(3-4) 


The  unknown  parameters  oP  the  desired  gamma  density  Punction 
are  estimated  using  the  method  oP  moments  [ReP.  33:p.  250] 
using  the  Pol  lowing  relations; 

_22  2 
r  =  est(p)  /est(<r„  )  and  X  =  est  (p ) /est  (a„  ).  (3-5) 

Using  this  approach,  the  gamma  density  Punction  used  to 

approximate  the  inter arrival  time  histogram  is  also  plotted 

on  Fig.  10.  Figs.  11  and  12  are  based  on  the  same  data  as 

Fig.  10.  However,  they  are  separated  according  to  diPPerent 

phases.  It  can  be  seen  Prom  the  r  and  X  parameters  in  Table 

III  that  there  is  a  signiPicant  diPPerence  between  phases  oP 


the  same  gap  noise  process  shown  in  the  histogram  even 
though  the  mean  value  oP  each  data  set  was  almost  the  same. 
The  data  in  Fig.  11  was  more  clustered  than  the  Fig.  12 
data.  This  results  in  a  higher  r  and  lower  X  parameter  Por 
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stogram  op  Interarrival  Times  (Case  1  -  Phase 


GAMMA  DENSITY  FUNCTION.  N=96 


Inter arrival  Times  (Case  1  -  Phase  B) 


the  same  mean  value.  Figs.  13  and  14  are  additional  sets  of 
gap  noise  data  that  were  observed  during  the  research. 

Table  III  summarizes  the  r  and  X  parameters  Fitted  to 


the  data  described  above. 


Case  1 


Case  1  Phase  A 


Case  1  Phase  B 


Case  2 


Case  3 


33.2 


63.4 


23.3 


44.3 


14.1 


37,000 
71 ,400 
25 , 900 
74,600 
19,200 


TABLE  III 


A  Final  model  parameter  will  be  called  T*  and  is  a  constant 
delay  or  oFFset  oF  every  second  pulse  group  that  exists  with 
respect  to  the  Fundamental  voltage  waveForm. 

In  summary  the  inputs  to  the  model  based  on  time-domain 
data  are; 

1.  An  estimate  oF  the  average  amplitude  oF  all 
observed  pulses. 

2.  An  estimate  oF  the  average  number  oF  pulses  per 
group  For  each  phase. 

3.  An  estimate  oF  the  mean  oF  the  interarrival  time 
oF  the  pulses  For  each  phase. 

4.  An  estimate  oF  the  variance  oF  the  interarrival 
time  between  pulses  For  each  phase. 

5.  An  estimate  oF  the  delay  in  the  start  time  oF  the 
pulse  groups  oF  one  phase  with  respect  to  the 
Fundamental  voltage  waveForm. 


DENSITY  FUNCTION 


The  average  amplitude  oF  pulses  From  diFFerent  gap  noise 
sources  is  highly  variable,  ranging  From  the  instrumentation 
noise  Floor  to  higher  than  any  observable  signal  in  the  HF 
band.  The  average  number  oF  pulses  per  group  has  been 
observed  to  vary  From  1  to  greater  than  20.  The  mean  oF  the 
interarrival  times  ranges  From  .1  to  1  ms  and  the 
coeFFicient  oF  variation  (standard  deviation  divided  by  the 
mean)  oF  the  interarrival  times  has  ranged  From  approx¬ 
imately  .1  to  .3  For  the  gap  noise  processes  we  have 
observed . 

D.  MODEL  DEVELOPMENT 

The  principal  goal  oF  this  chapter  is  to  describe  a 
noise  model  that,  with  suitable  choice  oF  parameters,  pro¬ 
vides  an  adequate  simulation  oF  the  actual  physical  noise 
mechanism  oF  gap  noise.  In  this  section  the  time-domain 
data  and  parameters  derived  in  the  previous  section  will  be 
integrated  into  the  Filtered  impulse  model  developed  in 
Chapter  II.  Recalling  From  Chapter  II,  a  Filtered  impulse 
noise  process  is  given  in  complex  envelope  notation  by 

N(t)  -j0, 

n(t)  =  g <t )  +  Z  a, e  h(t-t,  )  (3-6) 
i  =*  1 

where  p(a)  is  the  amplitude  distribution  oF  the  impulses, 
N(t)  is  a  unit  counting  process  whose  statistics  determine 
the  impulse  arrival  times,  0  is  a  random  variable  uniFormly 


distributed  between  0  and  2n  and  h(t)  is  the  impulse 


|  response  of  the  receiver  -filter. 

i 

'  In  order  to  physically  justify  a  filtered  impulse  model, 

i 


|  the  impulse  duration  must  be  small  compared  to  the  inverse 

i 

I  bandwidth  of  the  receiver  filter.  This  condition  is  easily 

met  in  the  case  of  gap  discharges.  Laboratory  analysis  of 
temporal  character  1 st i cs  of  gap  discharges  for  various 
|  geometries  shows  that  the  impulse  durations  range  from  10  to 

I 

100 's  of  nano-seconds  CRef.  341.  Therefore,  for  f i Iter 

bandwidths  up  to  1  MHz  the  output  noise  process  will  only  be 

a  function  of  the  incident  time  of  the  impulse  and  the 

filter  response,  not  the  waveform  of  the  impulse. 

The  specification  of  the  statistics  of  N(t),  the  count- 
t 


ing  process  that  drives  the  model  is  then  crucial  to 
obtaining  an  accurate  representation  of  the  physical  noise 
process.  Based  on  the  electrical  character i st l cs  of  a  gap 
discharge  and  the  observed  time-domain  characteristics ,  the 
following  assumptions  are  made; 


1.  A  primary  series  of  event  times  separated  by  an 
interval  T0 /2  exists.  T0  is  the  fundamental 
period  of  the  power  line  voltage  waveform.  The 
distribution  of  the  time  to  the  first  primary 
event  is  uniform  over  0  to  T0  . 


2.  A  subsidiary  process  commences  at  the  primary 
event  times  i.e.  T0 /2  ,  T0  ,  3T0 /2  ...  The 

subsidiary  process  is  a  renewal  process  that 
continues  for  N1  or  N2  renewals  alternating 
between  N1  and  N2  renewals  at  successive  primary 
po l nts . 


3.  The  interarrival  -times  to  -the  first  and  subsequent 
points  of  the  subsidiary  renewal  process  for  each 
phase  are  all  independent  and  identically 
distributed  (110)  within  a  pulse.  Each  phase  has 
its  own  set  of  parameters. 

4.  The  weighting  or  amplitude  distribution  of  all 
impulses  in  the  subsidiary  process  is  a  constant. 

5.  A  constant  offset  or  delay  designated  by  T»  may 
exist  at  every  second  primary  renewal  point. 

The  first  assumption  is  supported  by  the  periodicity  in 
the  data  that  is  related  to  the  fundamental  frequency  of  the 
power  line.  The  large  majority  of  observed  gap  noise  sources 
had  pulse  groups  on  both  the  positive  and  negative  phases  of 
the  fundamental  waveform.  This  accounts  for  the  T0 /2  per¬ 
iodicity. 

The  second  assumption  addresses  the  differences  between 
the  sparking  phenomenon  on  the  positive  and  negative  phases 
of  the  fundamental .  It  is  frequently  observed  that  the 
ECN»««]  is  different  than  ECN«aa3  where  N* « •  and  N*  a «  are 
the  random  number  of  discharges  in  the  pulse  groups.  This 
can  be  due  to  two  effects;  an  asymmetrical  gap  geometry  and 
the  fundamental  physical  difference  between  the  sparking 
mechanisms  for  positive  to  ground  and  negative  to  ground 
sparks.  Two  deterministic  integer  constants;  N1  and  N2 


nearest  to  ECNpoa3  and  ECNHaa3  respectively  are  assumed  as 
the  approximation  to  the  variable  number  of  sparks  per  half 
cycle  of  the  fundamental  waveform. 
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The  third  assumption  is  based  on  the  fact  that  the 
"inception  of  gap  discharges  in  natural  air  and  the 
development  of  electron  avalanche  are  fundamentally 
probabilistic  processes  that  depend  on  atmospheric  pressure, 
humidity,  presence  of  natural  ions,  electrode  surface  and  so 
on"  CRef.  34].  In  view  of  the  above  statement  and 
considering  the  empirical  data,  the  justification  for  devel¬ 
oping  the  gap  discharge  as  a  probabilistic  process  is  well 
founded.  The  assumption  of  independence  from  discharge  to 
discharge  is  not  as  well  justified.  Effects  such  as 
electrode  heating  after  the  initial  discharge  in  a  cycle 
could  act  to  make  the  average  interarrival  time  vary  from 
pulse  to  pulse  within  a  pulse  group.  However,  to  develop  a 
tractable  model  the  assumption  of  independence  between 
arrival  times  is  made.  The  gamma  density  function  deter¬ 
mined  by  the  estimate  of  the  mean  and  variance  of  the 
interarrival  times  for  each  phase  is  used  to  define  the  PDF 
of  the  interarrival  times. 

The  fourth  assumption  concerns  the  amplitude  of  the 
impulses.  In  most  cases  the  amplitude  of  the  impulses 
within  a  pulse  was  nearly  constant.  In  some  cases  there  was 
a  difference  in  the  average  amplitude  from  negative  phase  to 
positive  phase,  however  to  keep  the  model  tractable  the 
amplitude  of  all  impulses  is  assumed  constant. 

The  fifth  assumption  incorporates  an  effect  that  will 


explain  one  of  the  features  observed  in  the  power  spectral 


density  (PSD)  of  the  observed  process.  The  physical 
rationale  behind  this  assumption  is  that  there  may  exist  a 
difference  in  the  threshold  between  the  positive  and 
negative  phases,  thus  consistently  delaying  or  offsetting 
the  start  of  the  renewal  process  on  one  of  the  phases. 

The  above  assumptions  place  this  model  in  a  class  o-P 
processes  known  as  branching  processes  with  the  main  process 
being  a  degenerate  renewal  process  with  an  interarrival  PDF 
of 

p(t)  =  6 (t-Ta /2 )  (3-7) 

The  subsidiary  process  is  a  renewal  process  o-P  N1  or  N2 
points.  Consequently,  the  impulse  source  For  the  gap  noise 
model  is  completely  specified  by  the  following  parameters : 
a  -  Amplitude  of  the  impulses, 

N1,N2  -  Number  of  impulses  in  the  subsidiary  process. 

Pi ,p»  -  Mean  value  of  impulse  interarrival  time, 

Var(p, ) , Var (p« )  -  Variance  of  interarrival  times, 

T*  -  Offset  on  one  phase. 

To  complete  an  interference  scenario  N* B« r ,  /2  and  h(t)  must 
also  be  specified.  The  expression  for  the  branching  renewal 
process  is  then 

M(T)  N.  -j0.B 

n(t)  «  g(t)  +  Z  Z  a.  „  e  h ( t-t.  „ -T„  )  (3-8) 

m=l  n= 1 

Fig.  15a, b  shows  sample  realizations  of  the  envelope, 
E(t),  of  the  noise  process  defined  above. 


Figure  15.  Envelope  and  Counting  Function  of  n(t) 


The  Gaussian  term  was  set  to  0  and  h(t)  was  approximated  as 
a  decaying  exponential.  Note  that  in  comparison  to  an 
actual  noise  process  illustrated  in  Fig.  8  the  amplitude  of 
the  impulses  is  constant  with  a  Fixed  number  oF  impulses  on 
the  positive  and  negative  polarity  oF  the  line  voltage 
waveForm.  Fig.  15c  is  the  primary  counting  process,  M(T)  , 
driving  the  subsidiary  process  shown  in  Fig.  15b.  Fig.  15b 
includes  a  non-zero  T*  and  its  relationship  to  M(T)  is  shown 
in  Fig.  15c. 

Fig.  16  shows  a  realization  oF  the  envelope  oF  this 
process  For  N1  and  N2  *  3,  T  =30  ms  and  gamma  density 
Function  parameters  rl  and  r2  =  32  and  XI  and  \2  =  37,000. 
Random  interarrival  times  with  a  gamma  distribution  were 


approximated  by  summing  exponentially  distributed  random 
variables  [Ref.  353  .  The  impulse  amplitudes  were  normalized 
to  one  and  a  background  level  of  quadrature  Gaussian  noise 
at  N« B/2  =  .0001  was  added  which  simulates  either  receiver 
noise  or  high  density,  low  amplitude  impulsive  noise. 

E.  COMPARISON  OF  CALCULATED  AND  OBSERVED  DATA 

One  partial  description  of  a  noise  process  is  the 
amplitude  probability  distribution  (APD)  o-f  the  envelope. 
This  description  has  been  used  extensively  in  the  analysis 
of  noise  and  in  deriving  optimum  receivers.  For  the  N(t) 
specified  in  the  previous  section  this  calculation  would  be 
difficult.  A  second  noise  process  descriptor  is  the 
spectrum  of  the  noise  process.  This  descriptor  is  suited  to 
our  model  where  the  structure  of  the  noise  is  contained  in 
the  counting  process  that  drives  the  impulse  generation.  In 
addition  to  receiver  noise  performance  evaluation,  spectral 
analysis  can  also  be  used  for  noise  source  1  dent i f 1  cat  1  on 
and  isolation  [Ref.  363. 

In  order  to  determine  the  spectrum  of  the  envelope 
squared  of  the  branching  renewal  process  defined  by  Eqn. 
3-8,  the  spectrum  of  the  envelope  squared  of  a  non-branching 
renewal  process  is  determined  in  Appendix  D.  The  results 
obtained  in  Appendix  D  are  then  used  to  derive  the  spectrum 
of  the  the  branching  renewal  process  which  is  postulated  to 
model  gap  noise.  Starting  with  the  branching  renewal 


process  given  in  Eqn.  3-8,  and  assuming  g(t )  is  equal  to 


zero,  it  is  shown  in  Appendix  E  that  the  average  o-P  the 
Bartlett  estimate  o -P  the  PSD  is: 

2  5 

S(w,T)  =  |H«(u))|  /T  E  Term  (n)  (3-9) 

n=  1 

where  n  is  an  index  to  the  5  Terms  de-Pined  below.  These 
■Five  Terms  involved  in  the  summation  in  Eqn.  3-9  are  defined 
as  -Pol  lows : 

Term  1  -  The  inter-pulse  group  summations  between  pulses  in 

pulse  groups  with  N1  pulses 

Term  1  =  M2  SCu.Nl)  (3-9a> 

Term  2  -  The  inter-pulse  group  summations  between  pulses  in 

pulse  groups  with  N2  pulses 

Term  2  *  M2  5(u>,N2>  C3-9b> 

Term  3  -  The  intra-pulse  group  summations  between  pulses  in 

groups  with  N1  pulses 

cos (M2  2nu/u, )  -  1 

Term  3  *  <  -  -  M2  >  (Hui.Nl)  (3-9c) 

cos(2nu/u),  )  -  1 

Term  4  -  The  intra-pulse  group  summations  between  pulses  in 
groups  with  N2  pulses 

cos(M2  2n<u/ua  )  -  1 

Term  4  =  { -  -  M2  >  0(tu,N2)  (3-9d> 

cos  (2nu/iii.  )  -  1 

Term  5  -  The  intra-pulse  group  summations  between  pulses  in 

groups  with  N1  pulses  per  group  to  pulses  in  groups  with  N2 
pulses  and  vice  versa 
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cos(M2  2nui/ui.  )  -  1 

Term  5  =  2cos(  (ui(T.  /2-T0  )  )  -C - ><|/(0,  ,  N1  ,0,  ,  N2) 

cos  ( 2ndi/(d,  )  -  1 

<3-9e> 

where  0,  and  0,  are  the  characteristic  -functions  associated 
with  the  interarrival  times  for  each  phase.  The  parameters 
N1  and  N2  are  the  number  of  renewals  associated  with  each 
phase  and  3,  ft,  and  <4/  are  defined  below: 


3(0, N) 


N  -  N  |  0  ( ju>)  | 


|  1  -  0(jui)  | 


ft (0 , N)  = 


0(jUl)  (1  -  0(jUl)  ) 

2  Re-C -  >, 

2 

(1  -  0  (  jui)  )  (3-10) 


IvH-l  2 

|0(jai)  -  0(jw)  | 


|  1  -  0(jO)>  | 


(3-11 ) 


<4/(0,  ,  N 1  ,0.  ,N2)  = 


Nl+1  N2+1 

(0,  (j«)  -  0,  (jUl)  )  (0,  (  -  jui)  -  0.  (  —  j  ui )  ) 


(1  -  0,  (  j  ui )  )  (1-0,  (  -  j  ui )  ) 


(3-12) 


As  a  test  of  tho  model  the  PSD  of  the  gap  noise  process 
shown  in  Fig.  9  was  determined  using  a  spectrum  analyzer  and 


also  analytically  determined  using  only  the  model  parameters 
from  the  time-domain  data  and  the  methodology  of  the 
previous  section.  Fig.  17  is  a  computer  plot  of  the 


analytic  estimate  of  the  PSD  which  is  compared  with  the 
observed  PSD  computed  on  a  Wavetek  UA500A  spectrum  analyzer 
shown  in  Fig.  18. 


16 


Figure  18.  Observed  Estimate  o-f  Power  Spectral  Density  CPSD) 
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0.0  250.0  500.0  7SO.O  1000.0  1250.0  1500.0  1750.0  2000.0  2250.0 

FREQUENCY  -  HZ 

Figure  17.  Predicted  Pow#r  Spectral  Density  (PSD) 

Comparing  the  -two  PSD 'a,  it  is  seen  that  tha  analytical 
axpraasion  corractly  pradictad  all  oF  tha  signiFicant 
Faaturaa  aaan  in  tha  obaarvad  PSD. 

Tha  components  oF  tha  spectrum  dua  to  tha  diFFarant 
terma  ia  oF  intaraat  and  oFFera  aoma  inaight  into  the 
origina  oF  tha  Faaturaa  oF  tha  obaarvad  apectrum.  Tha  dom¬ 
inant  Feature  oF  tha  obaarvad  apectrum  are  tha  periodic  120 
Hz  apactral  linaa.  At  800  Hz  tha  apectral  linaa  change,  and 
appear  aa  odd  harmonica  oF  60  Hz,  however,  still  remaining 
at  120  Hz  interval  a.  Tarma  3-5  are  raaponaibla  For  tha 
thaaa  apactral  line  Faaturaa.  Each  oF  thaaa  tarma  is  a  comb 
Function,  determined  by  tha  Fundamental  Frequency  and 
observation  interval,  multiplied  by  an  envelope  determined 


by  the  inter arrival  time  distribution.  Figs.  19  and  2 0  are 
Terms  3  and  4  .  respectively  and  clearly  show  the  phase  to 
phase  diFFerence.  The  more  peaked  behavior  of  Fig.  19  at 
1100  Hz  is  due  directly  to  the  more  clustered  interarrival 
time  behavior  on  one  phase.  This  eFFect  Mas  observed  in  the 
time  domain  in  Fig.  11.  When  the  terms  shown  in  Figs.  19 
and  20  are  added  together,  they  Form  60  Hz  harmonic  spectral 
lines  and  when  Term  5,  shown  in  Fig.  21,  is  added,  positive 
re i nForcement  occurs  at  120  Hz  and  negative  at  60  Hz 
harmonics.  This  produces  the  observed  120  Hz  harmonic 
spectral  lines.  When  the  cancellation  is  not  complete  due 
to  diFFerences  in  the  character istic  Function . From  phase  to 


O 


FREQUENCY  -  HZ 

Figure  21.  PSD  Term  due  to  Phase  A  to  Phase  B  Intra-Pulse  Group 
Correlations 


phase,  small  60  Hz  components  can  be  present.  This  eFFect 


is  seen  in  both  the  observed  and  predicted  spectra. 

The  observed  spectra  changeover  from  120  Hz  harmonics  to 
odd  60  Hz  harmonics  at  800  Hz  is  also  correctly  predicted  by 
the  model.  This  eFFect  is  due  to  the  factor ,  T»  ,  appearing 


in  Term  5 

which  causes  the 

odd 

60 

Hz 

harmoni cs 

to 

be 

positively 

reinForced  and 

the 

120 

Hz 

harmonics 

to 

be 

canceled  out  above  800  Hz.  The  ■Pinal  e-P-Pect  shown  in  the 
model  and  seen  in  the  observed  spectra  is  the  presence  o-P 
continuous  spectral  components.  The  continuous  spectra  due 
to  Terms  1  and  2  is  shown  in  Figs.  22  and  23  and  also 
di-P-Pers  -Prom  phase  to  phase.  The  continuous  spectra  peaks 
markedly  at  1100  Hz  and  this  e-P-Pect  is  apparent  in  the 
observed  spectra  as  a  rise  in  the  noise  floor.  The  con¬ 
tinuous  spectral  term  arises  when  Term  1  is  added  to  Term  3 
and  Term  2  to  Term  4.  The  peak  value  at  1100  Hz  in  Terms  1 
and  2  is  greater  than  the  negative  values  in  Terms  3  and  4 
at  the  same  Frequency  and  when  added  together  causes  a 
continuous  rise  in  the  spectral  Floor.  The  reason  For  this 
is  that'  the  i  to  i+1  and  i-1  interarrival  times  within  a 
pulse  group  are  more  correlated  than  the  correspond i ng 
interarrival  times  between  diFFerent  pulse  groups.  The  more 
correlated  arrivals  appear  in  Terms  1  and  2  and  cause  the 
spectral  peak  to  be  larger. 
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Figure  22.  PSD  Term  due  to  Phase  A  Inter-Pulse  Group 
Correlations 
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Figure  23.  PSD  Term  due  to  Phase  B  Inter-Pulse  Group 
Correlations 


F.  SUMMARY 


strength  of  air,  which  is  about  30  kV/cm  -For  gaps  over  1  cm 


in  length,  prior  to  the  spark  breakdown  voltage,  then  corona 
will  exist.  Sharp  points  and  small  radius  wires  are  avoided 
whenever  possible  on  power  lines,  and  corona  generally  does 
not  become  a  significant  source  of  radio  interference  until 
voltages  exceed  about  70  kV  CRef.  28:p.  841. 

As  the  voltage  on  an  AC  line  is  increased,  corona  will 
appear  and  be  associated  with  either  the  positive,  negative 
or  possibly  both  half-cycles  of  the  fundamental  waveform. 
For  an  aged,  clean  transmission  line  it  appears  that  the 
threshold  voltage  for  negative  corona  is  less  than  that  of 
positive  corona.  If  the  transmission  line  has  any  unusual 
characteristics,  the  threshold  voltages  of  positive  and 
negative  corona  become  difficult  to  categorize  CRef.  371. 
In  particular,  it  has  been  observed  that  positive  streamer 
corona,  which  is  more  disruptive  of  communications  fre¬ 
quencies  than  negative  corona,  predominates  on  newly 
installed  lines  and  during  periods  of  precipitation.  For 
modeling  purposes,  the  important  character istics  of  corona 
noise  are  not  a  precise  description  of  the  physical 
mechanisms  but  the  relationship  between  the  observed  inter¬ 
ference  waveform  and  the  line  voltage. 

The  procedure  most  commonly  used  for  the  prediction  of 
interference  due  to  corona  noise  gives  the  noise  power  as  a 
function  of  the  characteristics  of  the  interference 
scenario.  The  random  properties  of  the  corona  noise  are 
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averaged  out  CReF.  383.  These  characteristics  would 
commonly  include  distance  From  the  line,  the  line  voltage 
and  a  weather  correction  term.  Each  oF  these  terms  is 
experimentally  veriFied  and  related  to  some  physical  process 
in  the  generation  oF  corona  noise.  While  this  procedure 
gives  an  estimate  oF  the  noise  power  in  a  specified 
bandwidth,  it  does  not  incorporate  the  probabilistic  nature 
or  time-domain  periodicity  observed  in  the  noise. 

The  approach  taken  in  this  chapter  will  continue  with 
the  Filtered  impulse  theme  oF  the  dissertation.  For  corona 
noise,  the  counting  process,  N(T)  ,  will  be  modeled  as  a 
periodically  modulated  stream  oF  impulses.  This  approach 
will  have  the  advantage  oF  being  accurately  able  to  simulate 
the  time-  and  Frequency-domai n  behavior  oF  corona  noise. 
However,  as  with  the  gap  noise  model  it  will  suFFer  From  the 
disadvantage  that  For  each  scenario,  a  set  oF  parameters 
must  be  exper imental ly  determined. 

C.  TIME  DOMAIN  CHARACTERISTICS 

The  time-domain  data  For  corona  noise  was  obtained  using 
the  instrumented  van  described  in  Appendix  A.  The  antenna 
was  a  3-meter  whip  mounted  on  the  van  (which  was  parked 
underneath  a  3-phase  500-kV  transmission  line).  The  weather 
conditions  at  the  time  oF  measurement  were  clear  and  sunny. 
Records  oF  the  corona  noise  were  obtained  using  the  envelope 
detected  output  oF  the  HP-141T  described  in  Appendix  A  and 


'V  'v.v 


;  v:  ■>  .  w; 777777  7:  ~  ~  \r  rrf: 


7  V 


■  y »  v* 


tXvTT> 


■  recording  -the  data  on  a  20-kHz  bandwidth,  85-dB  dynamic 

|  range  digital  audio  system.  The  recording  process  allowed 

I  subsequent  analysis  in  the  laboratory. 

A  37-ms  portion  oF  a  data  record  is  shown  in  Fig.  24. 
The  data  was  taken  at  a  300-kHz  center  Frequency  and  10-kHz 
intermediate  Frequency  (IF).  This  particular  picture  was 
|  obtained  by  playing  back  the  tape,  re-digitizing  the  data 

!  and  displaying  it  with  a  signal-analysis  software  package. 

>  One  oF  the  important  Features  oF  the  time-domain  record  in 
Fig.  24  is  the  periodicity.  The  Fundamental  Frequency  oF 

.  180  Hz  can  be  determined  From  a  periodic  rise  in  the  mean 

'  value  oF  the  noise  and  From  a  periodic  increase  in  high- 

amplitude  impulses.  The  Fundamental  noise  Frequency  oF  180 
j  Hz  on  a  three  phase  line  indicates  that  the  dominate  noise 

|  is  generated  on  one  halF-cycle  oF  the  Fundamental  60-Hz 

[  waveForm  on  each  oF  the  three  phases  oF  the  line.  This  is 

.  consistent  with  low-level  negative  corona  predicted  For  aged 

lines  in  clear  weather. 

The  periodogram  spectrum  oF  the  data  record  described 
,  above  is  shown  in  Fig.  25.  The  time-domain  observations  oF 

I 

»  the  180-Hz  periodicity  are  conFirmed  and  a  360-Hz  harmonic 

l 

>  oF  the  Fundamental  is  also  present.  At  higher  Frequencies, 

\  no  harmonics  are  observed  and  the  spectrum  becomes  white. 

! 

I  This  is  in  distinct  contrast  to  the  gap  noise  described  in 

I 

•  Chapter  III  where  the  spectrum  oF  gap  noise  has  harmonics 

1 

j  well  into  the  kHz  range.  The  periodicity  can  also  be  seen 
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Figure  25.  Corona  NQise  in  Frequency  Domain 


in  Fig.  26  (which  is  the  same  data  recorded  on  the  3-Axis 
display).  The  contribution  oF  each  phase  of  the  power  line 
to  the  total  corona  noise  process  can  be  clearly  seen  in  the 
lower  picture  in  Fig.  26. 

Another  important  Feature  oF  the  data  For  modeling 
purposes  is  that  the  receiver  response  to  an  individual 
impulse  overlaps  with  the  response  oF  preceding  impulses. 
Although  the  noise  is  still  impulsive,  the  simpliFying 
assumption  oF  non-overlapping  pulses  is  not  valid.  It  is 
seen  in  Fig.  24  that  the  amplitude  oF  the  impulses  is  random 
and  there  is  no  Fine-grain  interarrival  time  structure  other 
than  the  periodicity  discussed  above.  Both  oF  these 
observations  are  in  direct  contrast  to  gap  noise  where  the 
amplitude  From  impulse  to  impulse  is  almost  constant,  and 
there  is  a  deFinite  impulse  interarrival  time  structure. 
These  observations  provide  a  method  to  distinguish  between 
gap  noise  and  corona  noise. 

An  interesting  aspect  oF  noise  From  transmission  lines 
was  noted  while  taking  measurements  For  the  corona  model. 
It  was  diFFicult  to  obtain  a  recording  oF  power-line  noise 
in  the  vicinity  oF  high  voltage  transmission  lines  due  only 
to  corona  noise*.  In  many  cases,  a  gap  noise  component  was 
superimposed  on  the  corona  noise.  This  can  be  seen  by 
comparing  Figs.  26  and  27.  Fig.  27  is  a  view  oF  both  gap 
and  corona  noise.  The  corona  noise  amplitude  peaks  at  the  2 
pV  level,  and  the  amplitude  oF  the  gap-noise  impulses  is 


about  3.5  mV.  The  two  low«r  pictures  are  oP  identical  data 
with  the  threshold  varied  to  emphasize  each  noise  type.  The 
lower  picture  has  the  threshold  set  at  about  3  jiV 
eliminating  the  180-Hz  corona  noise  and  showing  the  gap 
noise  process  occurring  at  a  120-Hz  rate.  This  corresponds 
to  a  gap  noise  source  sparking  on  the  positive  and  negative 
halP-cycle  o-P  a  single  phase  o-P  a  3  phase  line.  When  the 
threshold  is  lowered  to  about  1  pV,  the  180-Hz  corona  noise 
dominates  as  seen  in  the  middle  picture  o-P  Fig.  27.  This 
phenomenon  is  due  to  corona  discharge  on  the  negative  halP- 
cycle  o-P  each  phase  o-P  a  3-phase  system.  Without  the  abil¬ 
ity  to  examine  the  time-  and  frequency-domain  behavior  o-P 
the  noise  using  a  scanning  analyzer  and  3-Axis  display,  it 
would  be  very  di-P-Picult  to  d i PPerent i ate  between  gap  and 
corona  noise  and  make  accurate  measurements  o-P  either. 

Fig.  28  illustrates  corona  noise  with  di-P-Perent  phase - 
to-phase  characteristics.  In  this  example  the  noise  power 
due  to  the  corona  is  unequal  Prom  phase  to  phase.  The  noise 
occurring  on  the  dominant  phase  is  approximately  twice  the 
peak  voltage  oP  the  phase  with  the  lowest  amplitude. 

D.  CORONA  NOISE  MODEL 

Based  on  the  above  observations  oP  the  demodulated 
envelope  oP  corona  noise,  a  speciPic  Piltered  impulse  model 
based  on  Eqn.  2-18  will  be  developed.  This  model  is  similar 
to  models  used  in  synthetic  hydrology  to  account  Por 


Figure  28. 


Corona  Noise  -  Uneven  Phase  to  Phase 


seasonal  variations  in  stream-flow  data  CRef.  391.  Assuming 
no  Gaussian  noise  component,  a  lowpass  complex  -filtered 
noise  process  can  be  written 

N<t>  -j0i 

n  <t )  »  E  a, e  & (t-t, ) .  (4-1) 

i-1 

Since  corona  noise  is  a  high  density  process,  the  envelope 
cannot  be  simply  de-fined  as  in  Eqn.  3-22  since  the  response 
o-f  the  -filter  to  adjacent  pulses,  will  overlap.  To  cir¬ 
cumvent  this  di-f-ficulty  it  will  be  assumed  that  the  envelope 
o-f  n(t)  can  be  modeled  by 

N  <  t ) 

e(t)  =  E  a,  g(t-t,>  (4-2) 

i  =  l 

where  both  e(t)  and  g(t)  are  real,  positive  and  g(t)  depends 
on  h<t)  and  the  process  intensity. 

The  assumptions  to  de-fine  the  model  are  as  -follows; 

1.  The  rate  parameter  -for  the  counting  function  N(t) 
is  periodically  modulated  at  -frequency  u>o  and  will 
be  approximated  by  a  -finite  number  o-f  terms  in  a 
Fourier  series. 

2.  The  amplitude  distribution  of  successive  pulses 
given  by  the  a,  's  are  independent  and  identically 
distributed  with  density  function  p(a) . 

3.  The  start  time  of  the  periodic  modulation  is 
uniformly  distributed  over  the  period  T0  where  T„ 

■  1  /w«  . 

Assumption  one  is  motivated  by  the  periodicity  observed  in 
the  data.  In  most  cases  the  fundamental  period,  u>o  will  be 
180  Hz.  If  the  noise  is  more  dominant  on  one  phase  of  the 
line,  as  shown  in  Fig.  28,  then  60  Hz  components  of  the 


<  s 


noist  will 


exist  and  -for  a  more  accurate  model  the 


■fundamental  period  would  be  60  Hz.  The  rate  parameter  is 
given  by 

Nt 

\(t)  *  £  a„  cos  (<i»o  nt )  (4-3) 

n=0 

where  Nt  is  the  order  at  which  the  Fourier  series  is 
truncated . 

The  second  assumption  is  arbitrary  in  terms  o-f  the 
distribution  Function  chosen  for  the  impulses.  It  is 
obvious  that  random  impulse  amplitudes  are  being  generated 
by  the  cumulative  effect  of  the  corona  noise  sources. 
However,  it  would  be  difficult  to  accurately  fit  a  density 
function  to  the  data  due  to  the  time  varying  nature  of  the 
process.  The  independence  of  the  amplitude  bursts  is  also 
conjecture  but  facilitates  the  development  of  a  tractable 
model.  The  Rayleigh  density  function  is 

2 

-a  / 2fl 

p( a)=(a/fl)e  a  2  0 

*  0  a  <  0  (4-4) 

and  was  picked  to  model  the  random  amplitudes.  It  is 
analytically  tractable  and  can  be  justified  by  assuming  the 
quadrature  components  of  the  impulse  amplitudes  are  in¬ 
dependent  Gaussian  events. 

The  final  assumption  is  made  because  in  most 
interference  scenarios,  the  power  line  noise  process  will 
not  be  synchronized  to  the  process  being  interfered  with. 
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Thus  the  phase  will  be  effectively  randomized  over  one 
•fundamental  period.  In  some  of  the  observations  of  corona 
noise  (see  Figs.  26-28)  this  assumption  was  clearly  violated 
and  the  instrumentation  was  intentionally  synchronized  to 
the  power  line  for  clarity  of  presentation. 

The  model  could  be  further  generalized  by  allowing  the 
amplitude  density  function  p(a)  to  be  a  periodic  function  of 
time,  p(a,t)  independent  of  the  rate  parameter .  This  added 
complexity  would  make  the  model  more  accurate;  however,  more 
parameters  would  have  to  be  estimated. 

The  above  assumptions  classify  this  process  as  a 
compound  non-homogeneous  filtered  Poisson  process.  A  con¬ 
venient  tool  to  study  these  processes  is  the  cumulant 
generating  function  CRef.  41  :p.  117] 

se(t> 

\y(s)  =  ln<ECe  ]>  (4-5) 
where  In  is  the  natural  logarithm.  The  joint  cumulant  gen¬ 
erating  function  is 

St  e  (ti  )  +  s»  e  ( t»  ) 

4/<St,St)  =  ln<ECe  ]>.  (4-6) 


For  the  noise  process  defined  by  the  first  two  assumptions, 
it  is  shown  in  Appendix  F  that 


i;/  (  s )  =  p  (  a ) 


sahtt-c ) 

X ( a )  C  e  -  1 Idada 


(  4-7  ) 


'  P  Si  ah(tt  -a)  +■  srah(t,  -a) 

4/(s,  ,s«)  =  p(a)  X(a)  C  e  -  13dada. 

-«  -«  ( 4-8 ) 


From  the  properties  of  the  cumulant  generating  Functions  it 


can  be  shown  that 


ECe(t) 3  =  4/' (0) , 
VarCe(t)3  =  4/''(0) 


(4-9) 


(4-10) 


Cov  (t,  ,  t»  ) 


2 

o  4/(0, 0) 


(4-11) 


os,  os* 


It  is  straight Forward  to  show  that 


ECe(t) 3 


EC  a  3  E  a„  cos  (<i>*  no  )  h  (t-a)  da  , 

J  n=  1 


(4-12) 


2  P  Nt  2 

VarCe(t)3  ■  ECa  3  E  a„  cos  (w«>  na  )  h  (t-a)da 

J  n=*  1 


(4-13) 


2  P  Nt 

CovCt,  ,ti  ]  *  ECa  ]  E  a„  cos  (wo  ncx  )  h (t,  -a)h(t«-a)da. 

n=  1 

-« 

(4-14) 

The  autocorrelation  Function  can  be  Found  in  terms  oF  the 
above  Functions  as 

RCt,  ,t«  )  =  Cov  (t,  ,  t,  )  +  ECe(t,  )  3ECe(t,  )3.  (4-15) 


«!.».*  «■  •  '  m  »  -  «  a  4  a" 


.-  .-  ->  ->  ■>  V  V  ■_■  ■>  •_*>  ■>  ’.a 


The  above  equation*,  prior  to  the  application  oF 
assumption  3,  de-Fine  a  wide  sense  cyclostationary  process 
CReF.  403  where 

ECeCt  *  mT® ) 3  *  ECe(t)3  (4-16) 

and 


R(tt  -*■  mT®  ,t«  +  mT®  )  »  R(t,  ,t«  )  .  (4-17) 
The  statistics  o-F  the  process  are  invariant  to  a  shi-Ft  of 
the  time  origin  by  integral  multiples  o-F  the  period,  T0  . 
This  character i zat ion  is  particularly  appropriate  when  the 
process  is  observed  in  synchronization  with  the  -Fundamental 
power  line  Frequency. 

When  assumption  3  is  applied,  the  process  start  times 
are  randomly  shi-Fted  over  one  period,  the  process  becomes 
wide  sense  stationary  and  the  mean  and  autocorrelation 
become 


and 


ECn(t) 3 


T® 

* 

ECn(t) 3  dt 

0 


(4-18) 


1 

R  (t  )  *  - 

T, 


The  power  spectral  density 


T® 

R ( t+T , t ) dt . 

0 

(PSD)  oF  the  process 


i  s 


(4-19) 
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R  (t  )  e 


d*T  . 


(4-20) 


The  •following  algorithm  is  an  example  of  how  the  noise 


process  specified  above  can  be  digitally  generated.  For  the 
algorithm  specified  below;  N0  is  the  number  of  samples  in 
one  period  To*  <n,  >  is  the  real  sequence  to  be  generated,  i 
is  the  index  variable  for  the  sequence  and  j  and  k  are  dummy 
index  variables. 

1  .  i  »  0 

2.  n< i )  »  0 

3.  Generate  a  Rayleigh  random  variable,  R,  with 
parameter  A  and  set  n(i>  =  n(i)  +  R 

Nt 

4.  Set  intensity  X  ■  E  a,cos(n>t  ni/N0  ) 

n=0 

5.  Generate  a  Poisson  random  variable,  P,  and  set  j  =* 
trunc(P/X  +  .5) 

6.  Set  n  <  k )  ■  0  for  k  »  i  +•  1  to  j 

7.  i  -  j 

8.  Go  to  Step  3 

The  above  algorithm  will  continuously  generate  a  modulated 
Poisson  impulse  process  that  may  be  digitally  filtered 
concurrent  with  the  impulse  process  generation  or  subsequent 
to  the  generation  of  a  complete  time  record. 

Two  additional  points  should  be  noted.  A  number  of 
impulse  amplitudes  could  stack  up  on  each  other  if  the 
truncation  operation  returns  a  zero  value.  This  corresponds 
to  the  physical  case  where  multiple  impulses  occur  and  are 
unresolved  within  one  sampling  period.  Before  a  record  of 
noise  samples  generated  by  the  above  process  is  used  a 
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truncated  uniform  random  variable  between  0  and  N0,  the  sampled 

period,  should  be  subtracted  -From  i.  This  ensures  -that  the 
start  time  is  uni -Form  ly  distributed  over  one  period. 


E.  PARAMETER  ESTIMATION 

Having  specified  the  model,  the  task  o-F  determining  the 
parameters  -For  a  given  interference  scenario  remains.  The 
parameters  to  be  determined  are  the  <a*>,  the  Fourier  Series 
coe-F-Ficients  o-F  the  periodically  varying  rate  Function  and 
the  fl  parameter  in  the  probability  density  Function  oF 
impulse  amplitudes.  As  an  example,  the  parameters  For  the 
noise  process  described  in  the  time  domain  observations 
Section  will  be  determined. 

One  method  oF  estimating  these  ‘ parameters  is  by 
considering  the  predicted  mean  and  variance  oF  the  process 
as  given  by  Eqns.  4-12  and  4-13  and  matching  these 
parameters  to  the  data.  In  order  to  do  this  analysis,  the 
corona  source  has  to  be  observed  as  a  cyclostat ionary 
process.  This  requires  synchroni z i ng  the  data  collection 
instrumentation  to  the  line  Frequency  so  that  the  same  point 
in  each  cycle  can  be  examined.  This  is  diFFicult  to  do  For 
Field  measurements. 

Another  option  For  determining  the  model  parameters  is 
the  predicted  stationary  mean,  variance  and  PSD  or  auto¬ 
correlation  Function.  The  PSD  as  deFined  in  Eqn.  4-20  is 
stationary  so  no  line  synchroni zat ion  is  required  and 


>' 

l1 


instrumentation  is  readily  available  to  make  a  periodogram 


estimate  oF  -the  PSD. 

For  the  measurements  taken  in  this  research,  the  lowpass 
impulse  response  of  the  Filter  Mas  approximately 

2 

-at 

h(t)  =  (a/n  )  e  t  2  0 

*  0  t  <  0.  (4-21) 


With  the  above  deFinition  the  Fi Iter  has  unity  gain  at  u 
equal  to  0.  The  eFFective  Filter  response  g(t)  was  approx¬ 
imated  by  h(t)  . 

In  Appendix  G  it  is  shown  that 

2  2 
-u«  n 

1/2  Nt  4a 

ECe (t)  1  *  (an/2)  C  e  aBcos(tu,  nt)  (4-22) 

n=0 

1/2 

where  (nfl/2)  is  the  First  moment  oF  p(a) ,  the  probability 

density  Function  oF  the  impulse  amplitudes.  The  stationary 
mean  is 

1/2 

ECe(t)3  *  (fln/2)  a,.  (4-23) 

The  covariance  oF  e(t)  is  shown  to  be 

2 

-a (t,  -t, )  /2 

Cov(t, ,t« )  ■  (a/n)2fle 


2 


'  Nt  -2au 

E  a,  cos  ( u,  n  (u+  ( t,  -t-t*  ) /2)  e  du  . 
n*0 


1 1,  -t,  | 


(4-24) 


Using  Eqns. 


4-15  and  4-19  and  substituting  -the  above 


results,  the  stationary  autocorre 1  at  1  on  function  is 

a 

—  err  /  2  1/2  2 

R(t)  =215  Ce  a<>  (a/8n)  +■  (1/4)  a0 


2  2 
-oio  n 


Nt 

•+•  ( 1  /  8 )  E  e 
n=  1 


4a  2 

a0  cos  ( uio  n-T  )  1  . 

(4-25) 


Taking  the  Fourier  transform,  the  resulting  PSD  is 


1/2  -ai  /2a  2 

S  ( oi )  =  215  C  ( 2 )  a«  e  (n  /2)a,  4(u) 


2  2 
-0io  n 


Nt 

(it  /8)  E  e 
n=  1 


4a  2 

a,  A  ( oi 


nuj0  )  . 
(4-26) 


The  terms  in  the  PSD  in  order  are;  (1)  a  white  noise  term 
multiplied  by  the  -Frequency  response  of  the  Gaussian  Filter, 
(2)  an  impulse  at  oi  =  0  From  the  mean  value  oF  the  process 
and  (3)  periodic  components  at  harmonics  oF  180  Hz  due  to 
the  modulation  oF  the  driving  Function. 

To  estimate  the  model  parameters,  the  sample  mean  and 
mean  square  estimates  were  used  in  conjunction  with  Eqns.  4- 
23  and  4-24  to  solve  For  a0  and  fl.  The  periodogram  estimate 
oF  the  spectrum  was  used  to  determine  the  magnitude  oF  the 


a,  and  a*  co-eFFi  c  l  ents  . 


The  Fourier  series  was  truncated 
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Figure  29.  Synthetic  Coron* 

estimating  the  model  parameters  was  developed.  The  ad¬ 
vantage  of  -this  modal  over  previous  work  is  -that  it  allows 
an  accurate  representation  o-f  the  time-domain  behavior  o-f 
corona .  Tha  disadvantage  o-f  this  modal  is  that  it  is 
scenario  dependent  and  for  each  i ntarferenca  situation  a  new 
set  o-f  model  parameters  must  be  specified.  One  limitation 
was  introduced  into  the  model  by  the  assumption  leading  to 
Eqn.  4-2.  The  lowpass  equivalent  filter  is  a  function  o-f 
the  bandpass  filter  and  the  process  intensity,  and 
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V. 


ENERGY  DETECTION  RECEIVERS 

A.  INTRODUCTION 

Energy  detection  receivers  in  the  high  frequency  CriF) 
radio  band  operate  in  an  impulsive  radio  noise  environment. 
To  predict  and  improve  the  performance  of  energy  detection 
receivers,  accurate  models  oF  the  signal  and  noise  envir¬ 
onment  must  be  used.  The  noise  models  developed  in  earlier 
Chapters  will  be  used  to  simulate  atmospheric  and  man-made 
radio  noise  environments.  Two  types  oF  energy  detection  re¬ 
ceivers  will  be  considered:  a  Fast  Fourier  transForm  (FFT; 
processing  system  and  a  compressive  receiver.  To  improve 
the  performance  oF  these  receivers  the  concept  of  the 
locally  optimum  receiver  is  introduced  to  suggest  a  robust 
post-detection  processor . 

B.  ENERGY  DETECTION 

The  energy  detector  receiver  measures  the  enerqy  in  3 
signal  over  a  specific  time  interval.  For  this  chapter  we 
will  consider  the  signal  to  be  a  modulated  sinusoid,  of  un¬ 
known  Frequency  at  the  receiver.  The  signal  is  then 
represented  as 

s(t)  =  Acos(u).t  9(t>>  (5-1) 
where  A  and  <ds  are  random  variables  and  8(t)  is  a  slowly 
varying  Function  so  that  the  signal  power  is  effectively 
confined  to  a  narrow  bandwidth  B. 


\  \  v  v  %  v'-  \  \  •  %  v.  \  v  wr  v\%  .v.r* -jv  -  y.%  v  v  v .viv.v.v.vkvkv.^ .% .»  +  *  r* k"* 


Defining  -the  problem  as  a  hypothesis  test 
H„  :  x(t)  =  n(t)  noise  alone 


H, :  x(t)  =  n(t)  +  s(t)  signal  plus  noise. 


(  5  -  3  ) 


The  well  known  energy  detection  receiver  [Ref.  42]  uses  f £ c  . 
the  received  energy  as  the  test  statistic  such  that 


T 

f  2 

Under  H«  :  t* »  =  n  (t)dt 


Under  H, :  t« a  =  Cn(t)  +  s(t)]  dt . 


(5-2 


For  a  sampled  bandpass  process  the  test  statistic  is 

N  2  2 

H0  s  t®  o  —  £  C  nc  t  n£  ,  ] 

1*1 


n  2  ; 

H,  :  t«t  =  [  Cnc,  +  sc,  ]  +  [n,  ,  +  s*.  ] 

1  =  1 


xc  .  =  nCl  +  sc 


or  xc,  =  nc  , 


and  likewise  for  the  quadrature  component. 

Fig.  30  shows  the  block  diagram  of  a  receiver  that  cal¬ 
culates  this  statistic.  This  receiver  will  be  called  a 
square  and  sum  receiver.  The  square  and  sum  receiver  is  the 
optimum  receiver  for  detecting  random  phase  signals  in 
Gaussian  noise  at  small  signal  levels.  At  large  signal  to 
noise  ratios  (SNR)  this  detector  is  very  close  to  optimum 


Figure  3®.  Square  and  Sum  Receiver 


[Ref.  9:p.  3703.  Whan  N  *  8T'  and  the  signal  is  sampled  at 
the  Nyquist  rate,  the  test:  statistic  is  the  energy  in  tne 
signal  in  a  period  T. 

When  it  is  desired  to  implement  the  receiver  shown  in 
Fig.  30  over  a  wide  frequency  band  relative  to  the  target 
signal  modulation  bandwidth  B,  a  number  of  options  may  be 
pursued.  The  channelized  receiver  approach  in  which  an 
independent  receiver  operates  at  each  frequency  will  often 
be  too  expensive.  Two  other  possible  methods  for  generating 
the  test  statistic  are;  <1)  a  compressive  receiver  as  a 
spectrum  analyzer  [Ref.  433  and  (2)  a  fast  Fourier  transform 
(FFT)  signal  processor. 


It  has  been  shown  that  an  FFT  signal  processor  exactly 
implements  the  test  statistic  given  in  Eqn.  5-3  For  the 
assumed  signal  in  a  single  FFT  Filter  CReF.  44].  Whether 
successive  samples  will  be  independent  or  not  depends  on  the 
number  oF  transForms  per  second  relative  to  the  signal 
bandwidth  B.  IF  the  sample  interval  is  greater  than  1/B 
then  successive  samples  will  be  independent. 

In  the  compressive  receiver,  shown  in  Fig.  31,  a  Fast 
sweeping  local  oscillator  linearly  scans  the  band  oF 
interest.  Any  narrow  bandwidth  signal  in  the  band  will 
appear  at  the  mixer  output  as  a  linear,  Frequency  modulated 
(FM)  ,  chirp  signal.  This  chirp  signal  is  passed  through  a 
weighting  Filter  to  reduce  sidelobes  and  into  a  dispersive 
delay  line  (DDL)  which  is  matched  to  the  inverse  oF  the 
linear  FM  sweep.  The  output  oF  the  DDL  is  envelope  detected 
and,  when  referred  to  the  sweep  time,  it  provides  an 
estimate  oF  the  spectrum  oF  the  input  signal.  It  is  easily 
shown  that  For  a  sinusoid,  the  detector  output  is 
proportional  to  the  signal  magnitude.  Thus,  the  compressive 
receiver  approximates  the  square  root  oF  the  ith  value  oF 
the  test  statistic  in  Eqn.  5-3.  The  samples  will  be 
independent  iF  the  scan  revisit  time  is  greater  than  1/B. 

C.  SQUARE  AND  SUM  RECEIVER 

In  this  Section  the  energy  detection  perFormance  oF  the 
square  and  sum  receiver,  shown  in  Fig.  30,  will  be  examined 


•the  CDFs  oF  the  signal  plus  noise  densities  were  not  Found. 
The  P0  is  plotted  as  a  Function  oF  SNR  For  the  same 
parametric  conditions  as  used  in  Fig.  32.  As  anticipated 
the  Pd  decreases  dramatically  as  the  noise  becomes  more 
impulsive.  In  the  case  where  VD  equals  12,  the  Pd  decreases 
From  greater  than  99  percent  to  1  percent  at  a  5  dB  SNR  For 
the  same  PP  A  . 


A  receiver  that  is  optimum  in  Gaussian  noise  suFFers  a 
major  degradation  when  the  noise  statistics  become  non- 
Gaussian.  One  possible  approach  to  improving  the  per- 
Formance  is  to  design  a  receiver  that  is  based  on  the 
statistical  parameters  oF  the  noise.  A  parametric  approach 
will  be  discussed  in  the  next  section. 

D.  LOCALLY  OPTIMUM  RECEIVER 

The  concept  oF  a  locally  optimum  receiver  provides  a 
general  methodology  For  nonlinear  receiver  design  in  the 
presence  oF  non-Gaussian  noise.  The  approach  is  to 
determine  a  receiver  structure  that  is  optimized  For  the 
diFFicult  small  signal  case  and  then  examine,  and  lF 
necessary,  mod l Fy  the  structure  For  the  large  signal  case. 
In  this  Section  the  locally  optimum  receiver  For  energy 
detection  in  the  presence  oF  Hall  type  impulsive  noise  will 
be  determined  under  the  assumption  oF  weak  signal  levels. 
Previous  work  has  considered  known  signal  detection  in 
Cauchy  noise  which  is  related  to  one  oF  the  Hall  models 
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where  p„ (x|H0 )  and  pt  (x| H, )  are  the  probability  density 
functions  -for  the  observation  vector  x  under  H„  and  H, 
respectively.  Since  the  noise  samples  are  independent  and 
identically  distributed,  pa(x|H.)  and  Pi<x|Ht)  can  be 
written  as  products  of  the  univariate  noise  densities: 


p„  ( x )  =  11  pi  ( x,  > 

i  =  1 


p,  (x)  =  n  pi  (x,  -  Asi  )  . 

1  =  1 


(5-5) 


(5-6) 


The  test  can  then  be  written 


A  ( x ;  A )  = 


n  pi  ( x i  -  Ast  ) 
i  =  1 


n  pi  (x,  ) 

1  =  1 


>  T  decide  H„ ,  or 


<  T  decide  H, 


C5-?) 


In  deriving  the  locally  optimum  receiver  several 
approaches  have  been  used.  Rudnick  [Ref.  473  derived  a  form 
o-P  the  locally  optimum  receiver  by  expanding  the  likelihood 
ratio  in  a  Taylor  series  about  x  and  truncating  after  the 
second  order  terms.  Capon  CRef.  483  showed  that  the  locally 
optimum  detector  maximizes  the  slope  o-P  probability  of 
detection  versus  signal  strength  function  at  the  origin  and 
employs  a  statistic  based  on  the  following  test 


■oA  (  x  ;  A  ) 


tio  (x)  = 


(5-8) 


A  =>  0, 


m 
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In  Appendix  H  it  is  shown  that  the  locally  optimum 
receiver  For  the  detection  oF  a  random  sinusoid  is 


t'm  (r)  =  r  g  (r,  ) 
1  =  1 


(5-9) 


where  g(-)  is  the  zero-memory  nonlinearity  (ZMNL)  dePined 


be  1  ow 


g  ( r )  = 


F'  '  (r)  F'  (r) 


rF(r) 


(5-10) 


and  r  is  a  vector  o-P  envelope  samples.  The  function,  P(r), 
For  the  Hall  model  o-P  atmospheric  noise  discussed  in  Chapter 
1 1  is 


-P(r)  = 


(5-1  1  ) 


2.  2  (Q-D/2 

(r  +  y  ) 


where  K  is  a  constant.  From  Eqn.  5-10  the  ZMNL  of  the  test 
statistic  is 


g  ( r )  = 


(0+1 ) (0+3) r 

2  2  2 
( r  +  y  ) 


2(0+1 ) 

2  2 
( r  +  y  ) 


(5-12) 


The  receiver  that  implements  this  test  is  shown  in  block 
diagram  Form  in  Fig.  34.  The  ZMNL ' s  For  the  unit  root-mean- 
square  Hall  models  (see  Table  II)  parameter i zed  by  V0  are 
shown  in  Figs.  35  and  36.  The  ZMNL  receiver  character i st i c 


For  values  oF  V6  equal  2.1,  3,  6,  10  and  14  are  plotted. 

Recalling  that  the  nonlinearity  is  optimized  For  van¬ 


ishingly  small  signals  the  ZMNL  character i st l c  may  be 


Figure  34.  Locally  Optimum  Receiver  for  Enargy  Detection 

intuitively  analyzed.  Small  values  of  the  sampled  envelope 
will  be  given  a  negative  weighting,  intermediate  values  are 
emphasized  and  large  values  of  the  sampled  envelope  are 
effectively  nulled.  *  As  the  noise  becomes  more  impulsive, 
which  corresponds  to  a  larger  value  of  V„  and  smaller  9,  the 
above  features  become  more  pronounced.  This  is  seen  in  Fig. 
35  where  the  2MNL  has  large  negative  values  that  extend  off 
the  plot  for  very  small  envelope  values.  The  maximum  value 
of  the  2MNL  occurs  at 

1/2 

x  =  y< (9+5) / (0  +  1 ) )  c  5- i 3  > 

and  the  zero  crossing  at 

1/2 

x  =  y (2/ ( 0+1 ) )  .  (5-14) 

Both  of  the  points  are  proportional  to  the  mean  value  of  the 


Fig.  37  is  the  Pr  *  versus  threshold  For  the  locally 
optimum  receiver  shown  in  Fig.  34.  The  number  N  of  samples 
is  equal  to  10.  The  Hall  noise  samples  are  parameter i zed 
with  the  same  V0  values  used  in  Fig.  32.  The  threshold  data 
obtained  From  this  plot  was  used  to  generate  the  P&  curves 
shown  in  Fig.  38.  The  only  complete  data  set  was  obtained 
For  the  Hal  12  case  with  equal  to  2.1.  Some  data  was 
obtained  For  Ve  equal  to  5.  The  peaks  in  Po  For  the  more 
impulsive  noise  occur  at  smaller  SNR  levels,  not  shown  on 
the  plot.  The  data  points  are  plotted  on  the  graph  since 
the  dashed  curves  only  approximate  the  true  values. 

It  can  be  seen  that  the  locally  optimum  receiver  does 
improve  the  P0  (For  the  same  PF * )  in  the  small  signal  case. 
In  comparing  the  values  with  Fig.  33  For  Ve  equal  to  2.1  the 
locally  optimum  receiver  has  a  perFormance  better  than  the 
square  and  sum  receiver  From  -10  dB  to  7  dB  SNR.  Above  7  dB 
SNR,  the  P o  For  the  locally  optimum  receiver  drops  to  0. 
This  is  due  to  the  Fact  that  the  receiver  is  optimized  For 
small  signals.  The  Hal  12  ZMNL  in  Fig.  36  shows  the  reason 
For  this.  It  peaks  at  approx imate ly  1.3  times  the  RMS  value 
oF  the  noise  and  approaches  0  as  the  input  goes  to  inFinity. 
This  behavior  oF  the  ZMNL  must  be  modiFied  to  obtain  a 
receiver  that  will  also  work  in  the  large  signal  case.  In 
the  next  Section  a  robust  receiver  structure  is  proposed 
that  does  this.  A  robust  receiver  is  d i FFerent l ated  From  a 


parametric  receiver 


in  that  the  robust  rece l ver  does  not 
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SIGNAL  TO  NOISE  RATIO  C.rS) 

Po  for  Locally  Optimum  Receiver  in  Hall  Noise  as 
Function  of  V0  For  P,  «  =  1  x  l<Dm* 


depend  on  the  parameters  of  the  noise  density  and  it  per¬ 
forms  well,  1  f  subopt  imal  ly ,  over  a  class  or  even  differing 
classes  of  noise  densities. 

E.  ADAPTIVE  LIMITER 

The  ZMNL  derived  in  the  previous  Section  is  a  relatively 
complex  -function.  Furthermore ,  it  does  not  work  for  large 
signal  levels.  In  this  section  we  consider  a  practical 
implementation  of  the  ZMNL  consisting  of  an  adaptive  limiter 
to  approximate  the  ZMNL  function  with  a  modification  to 
improve  large  signal  performance.  Although  the  locally 
optimum  receiver  was  derived  for  Hall  type  impulsive  noise, 
if  the  approx  1  mat  1  on  to  it  is  properly  designed,  it  should 
show  robust  performance  in  the  presence  of  many  impulsive 
noise  process. 

Examining  the  P0  curves  Of  the  locally  optimum  receiver 
for  the  Hal  12  model  shows  that  it  provides  increased 
detection  performance  over  the  square  and  sum  receiver  from 
-10  dB  to  7  dB  SNR.  Since  this  region  is  where  the  square 
and  sum  receiver  suffered  the  worst  performance  degradation 
in  impulsive  noise,  the  Hall2  ZMNL  was  chosen  to  be  modified 
to  obtain  a  robust  nonlinearity.  The  solid  curve  in  Fig.  39 
is  the  Hal  12  ZMNL.  To  correct  its  performance  at  large 
signal  levels,  the  ZMNL  was  first  altered,  as  shown  by  the 
dashed  line  in  Fig.  39,  by  extending  its  character i st i c 
straight  out  from  its  maximum  value.  Furthermore,  an 


equivalent  test  statistic  is  obtained  if  g'(x)  is  any  linear 
function  of  g(x)  CRef.  461  so 


g'(x)  =  ag(x)  +6  a  ^  0. 


(5-15) 


Figure  39.  Modification  to  ZMNL 

With  suitable  parameters  a  and  3,  the  modified  ZMNL  is  shown 
in  Fig.  40.  The  level,  a,  remains  to  be  determined. 

Recalling  Eqn.  5-13,  the  argument  of  the  ZMNL  char¬ 
acteristic  at  its  maximum  value  is  proportional  to  the  mean 
value  of  the  Hall  distribution  for  which  it  was  derived  (see 
Table  I).  For  the  Hall2  distribution  the  proportional lty 
constant  is  approximately  1.6.  Therefore ,  in  order  to 

preserve  the  improvement  in  PD  in  the  same  SNR  region  as  the 


Hal  12  ZMNL,  the  limiter  level  was  set  to  equal  1.6  times  the 


Figure  40.  Linearized  Modification  to  2MNL 


Note  that  the  ZMNL  must  be  adjusted.  This  information  could 
be  known  a  priori  or,  as  symbolically  shown  in  Fig.  41, 
determined  from  the  statistics  of  the  noise.  Both  cases 
will  be  examined  below. 

The  false  alarm  rate  performance  of  this  modified 
receiver  in  Hall  and  Gaussian  noise,  where  the  statistics 
are  known  a  priori,  is  shown  in  Fig.  42  for  N  equal  to  10. 
The  value  for  a  was  determined  from  the  noise  mean.  At  a 
Pf *  equal  to  1  x  10" * ,  the  P0  curves  of  this  receiver  are 
shown  in  Fig.  43.  The  threshold  was  determined  from  Fig. 
42.  The  Gaussian  noise  performance  of  this  receiver  is 
poorer  than  that  of  the  square  and  sum  receiver  but  only  by 
a  few  percent.  However,  there  is  a  dramatic  improvement  i  r, 
the  PD  for  the  impulsive  noise  cases.  It  is  interesting  to 
note  that  for  the  same  noise  power,  a  well  designed 
nonlinear  receiver  provides  better  detection  performance  in 
impulsive  noise  than  the  optimum  receiver  for  Gaussian  noise 
CRef.  491. 

By  plotting  the  empirically  determined  threshold  from 
Fig.  42  versus  the  mean  of  the  noise  distribution  it  is  seen 
that  the  relationship  is  almost  linear.  This  is  shown  in 
Fig.  44  for  two  false  alarm  rates.  Using  this  relationship 
the  ZMNL  and  threshold  can  be  adaptively  modified  to  provide 
nearly  constant  false  alarm  rate  (CFAR)  performance  under 
varying  noise  conditions.  This  is  a  desirable  charac¬ 
teristic  for  an  energy  detection  receiver. 


A  -Final  question  to  be  resolved  in  lmplementinq  the 
receiver  is  the  real  time  determination  oF  an  estimate  of 
the  mean  and  mean  square  values  oF  the  noise.  One  approach 
is  to  sample  the  noise  in  a  signal  -Free  adjacent  channel  and 
another  is  to  -Form  the  sample  mean  oF  the  noise  using  past 
sample  values  with  no  signal  present  CReF.  501.  We  will 
implement  the  latter  approach  and  test  the  adaptive  receiver 
on  actual  signal  and  noise  samples.  In  sequencing  the 
samples  two  general  approaches  may  be  used.  The  samples  may 
be  processed  in  blocks  with  no  overlap  which  is  the  case 
correspond  1  ng  to  +-he  hypothesis  tests  discussed  in  this 
Chapter.  The  other  option  is  sequentially  processing  each 
sample  with  rank  two  updates  oF  the  test  statistic,  the  2MNL 
and  the  threshold.  In  this  case  the  hypothesis  test 
Formulation  will  provide  a  lower  bound  on  P6 . 

The  sequential  sampling  approach  was  chosen  and  two  runs 
oF  the  adaptive  receiver  are  shown  in  Figs.  45  and  46.  The 
test  statistic,  shown  by  the  dotted  line,  was  Formed  From 
the  10  previous  samples  using  the  receiver  structure  in  Fig. 
41.  The  adaptive  threshold,  shown  by  the  dashed  line,  was 
Formed  From  the  sample  mean  and  mean  square  estimates  oF  the 
previous  20  though  50  samples.  The  sample  mean  is  shown  by 
the  heavy  solid  line  and  the  individual  envelope  data  points 
by  the  light  solid  line. 

This  receiver  structure  requires  H0  to  be  true  at  the 
start  oF  the  test  and  will  detect  the  transition  to  H, . 


Since  the  adaptive  threshold  assumes  noise  only  samples,  it 
will  be  corrupted  when  signal  plus  noise  samples  are  used  to 
•form  it.  This  will  occur  20  samples  after  H,  becomes  true. 

This  effect  can  be  seen  in  Fig.  45  where  an  actual 
signal  from  one  filter  of  a  compressive  receiver  is  plotted. 
The  noise  is  highly  impulsive  and  would  be  a  very  difficult 
detection  environment  for  a  receiver  optimized  for  Gaussian 
noise.  Samples  50  to  110  are  noise  only  and  samples  110  to 
220  are  signal  plus  noise.  The  test  statistic  exceeds  the 
threshold  at  sample  125  and  the  receiver  would  declare  H, 
true  at  this  point.  However,  due  to  the  fact  that  the  test 
statistic  and  the  threshold  are  both  corrupted  by  the  signal 
plus  noise  samples,  the  test  statistic  later  drops  below  the 
new  threshold. 

In  Fig.  46,  the  receiver  performance  in  a  simulated  gap 
noise  environment  is  analyzed.  The  signal  commenced  at  sam¬ 
ple  110  and  the  receiver  declared  the  signal  present  at 
sample  117.  Again  the  threshold  is  1 ater  corrupted  and  the 
test  statistic  drops  below  the  threshold. 

F.  SUMMARY 

This  Chapter  has  presented  a  systematic  development  of 
the  bandpass  energy  detection  problem  in  impulsive  noise. 
Starting  with  a  Hall  density  function  for  the  noise,  the 
locally  optimum  receiver  is  derived.  This  receiver  is 
modified  to  design  a  practical  receiver  with  robust 


performance  in  non-Gaussian  noise.  The  design  methodology 
clearly  delineates  the  role  of  the  test  statistic,  threshold 
and  nonlinearity.  This  allows  independent  evaluation  of 
each  of  these  features.  One  aspect  of  energy  detection 
receivers  not  considered  is  the  role  of  data  in  filters 
adjacent  to  the  f i Iter  being  processed.  In  an  impulsive 
noise  environment,  the  noise  statistics  of  adjacent  filters 
will  be  correlated  and  their  use  may  potentially  offer  a 
further  performance  gain. 


V  I .  CONCLUSION 


A.  RESULTS 

The  dissertation  has  investigated  -the  modeling  of 
atmospheric  and  man-made  radio  noise  and  applied  the  results 
to  the  bandpass  energy  detection  problem. 

A  generalized  means  of  describing  an  impulsive  noise 
process  in  terms  of  its  complex  envelope  was  derived. 
Previous  models  of  atmospheric  noise  were  surveyed  and  some 
extensions  to  the  Hall  model  -For  atmospheric  noise  were 
developed. 

Using  the  physical  character  1  st  1  cs  of  gap  noise,  a  nine 
parameter  model  was  constructed  that  allows  arbitrary  sour¬ 
ces  o-P  gap  noise  to  be  synthesized.  The  driving  impulses 
-For  the  gap  noise  were  postulated  to  be  points  of  a 
branching  renewal  process.  The  application  of  this  type  of 
probabilistic  model  to  power  line  noise  is  unique. 
Additional  results  include  a  new  proof  of  the  power  spectrum 
of  an  equilibrium  renewal  point  process  and  a  derivation  of 
the  power  spectrum  of  a  truncated  branching  renewal  point 
process.  The  postulated  model  was  shown  to  correctly 
predict  all  the  significant  features  of  the  observed 
spectrum . 

In  a  similar  fashion,  corona  noise  was  postulated  as  a 
filtered  impulse  process  with  a  periodically  modulated  rate 
function.  This  model  was  shown  to  correctly  duplicate  the 


de-Fine  the  range  oF  values  For  the  power  line  noise  model 
parameters  and  to  determine  the  eFFect  oF  power  line 
construction  practices  on  the  parameters.  The  suscept¬ 
ibility  oF  various  types  oF  communication  equipment  to 
simulated  gap  and  corona  noise  should  also  be  examined.  The 
corona  noise  model  needs  to  be  extended  so  that  it  can  be 
speciFied  and  estimated  using  a  bandpass  vice  lowpass  Filter 
impulse  response.  The  work  on  energy  detection  suggests 
many  opportunities  For  Further  research.  The  problem  oF 
estimating  the  unknown  noise  parameters  For  the  adaptive 
updates  and  the  inclusion  oF  adjacent  Fi Iter  samples  into  an 
algorithm  are  prime  examples. 


APPENDIX  A 


INSTRUMENTATION 


The  instrumentation  configuration  employed  -to  provide 
data  on  the  detailed  time-  and  -frequency-domain  properties 
of  high  frequency  (HF)  radio  noise  is  shown  in  Fig.  47.  The 
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Figure  47.  Simplified  Slock  Diagram  of  Measurement  System 

radio  frequency  (RF)  input  to  the  instrumentation  was  either 
a  fixed  HF  long  wire  antenna  or  a  whip  antenna  mounted  on  a 
mobile  van.  One  complete  measurement  system  is  trans¬ 
portable  and  may  be  moved  from  site  to  site.  The  other 
system  is  installed  in  a  mobile  van  for  noise  measurements 
directly  at  the  interference  source  CRef.  531. 

In  the  lower  half  of  Fig.  47,  a  Hewlett-Packard  141T 
Spectrum  Analyzer  is  used  as  a  scanning  receiver  -to  drive  a 
Develco  7200B  3-  Axis  Display.  The  spectrum  analyzer  can  be 


■tuned  "to  any  desired  -frequency  in  -the  HF  band.  Its  scan 
rate,  scan  width,  intermediate  -frequency  (IF)  bandwidth,  IF 
gain,  RF  attenuation,  and  other  controls  can  be  adjusted  to 
best  describe  the  noise  under  observation. 

An  alternate  and  complementary  part  o-f  the  system  is 
used  to  examine  the  narrowband  properties  of  HF  noise  and  is 
depicted  in  the  upper  half  of  Fig.  47.  An  HF  receiver  is 
used  as  an  ampl if 1 er/trans 1 ator  and  tuned  to  a  frequency 
where  the  noise  is  present.  The  demodulated  audio  output  is 
then  applied  to  a  Wavetek  (Nicolet)  UA500A  spectrum  analyzer 
which  subdivides  the  audio-output  spectrum  into  about  500 
frequency  segments.  When  the  full  audio-output  bandwidth  of 
5  kHz  is  examined,  the  UA500A  provides  a  frequency 
resolution  of  10  Hz  per  segment.  The  Wavetek  analyzer  can 
provide  individual  transforms  or  average  a  selected  number 
of  successive  transforms.  The  transforms  (either  individual 
or  averaged)  are  then  presented  on  the  3-axis  display  or  an 
osc ill oscope . 

The  3-axis  display  provides  a  continuous  moving  real¬ 
time  visual  representat i on  of  the  analog  output  from  the 
spectrum  analyzers  or  narrowband  receiver  audio  output.  The 
receiver  or  analyzer  output  analog  data  are.  digitized  in  the 
display  and  stored  in  a  semiconductor  memory.  The  data  in 
the  memory  are  formatted  and  shown  on  a  cathode  ray  tube 
(CRT)  in  a  convenient  f requency - amp  1 i tude - 1 1 me  (3-axis) 


■format .  The  3-axis  presentation  can  be  frozen  at  any  desired 


time  and  photographed  with  a  standard  oscilloscope  camera. 

Fig.  48  illustrates  the  procedure  used  to  make  the  3- 
axis  presentation.  The  analog  input  is  divided  into  512 


ANO/OR 

SCAN  TIME  AXIS 

Figure  48.  Diagram  of  Data  Format 

equally  spaced  data  points  (indicated  by  the  horizontal  dots 
of  Fig.  48).  The  signal  amplitude  at  each  data  point  is 
represented  by  an  8-bit  word.  When  a  scan  is  completed,  its 
data  are  also  stored  in  memory.  Line  1  in  the  view  moves  to 
line  2,  and  the  new  scan  appears  as  line  1.  Subsequent 
scans  move  earlier  data,  line  by  line,  upwards  along  the 
time  axis  to  create  a  rising  raster  type  display.  When  the 
memory  is  full  (64  scans),  each  new  scan  is  entered  into  the 
bottom  line,  and  the  oldest  data  at  the  top  line  is 
discarded.  The  resulting  animated  view  provides  a  visual 
picture  of  noise  and  signals  within  the  block  of  frequencies 
if  the  spectrum  analyzer  output  is  being  observed.  If  the 


receiver  audio  output  is  being  observed  "the  display  is  a 
stacked  series  of  consecutive  time  records. 

The  3-axis  display  has  a  number  of  controls  to  aid  the 
observer  in  analyzing  the  structure  o-f  data  presented. 
Among  these  controls  are:  Cl)  elevation  and  azimuth  geometry 
controls  to  vary  the  viewing  aspect,  (2)  an  amplitude 
compression  control,  C3)  a  threshold  control  to  vary 
background  levels,  (4)  time-axis  expansion  controls,  and  (5) 
a  stop-action  switch  to  freeze  the  data  in  memory  for 
detailed  observation,  for  observation  from  various  aspects, 
or  for  photographing.  These  controls  are  used  to  optimize 
the  presentation  of  desired  signal  detail. 

Accurate  frequency,  time  and  amplitude  calibrations  are 
maintained  so  that  the  resulting  3-axis  photos  can  be 
manually  scaled  for  precise  signal  detail.  The  digitized 
data  at  the  display  input  are  available  at  a  digital 
interface  connector  for  external  digital  recording  or 
processing.  Received  data  can  be  processed  by  computer  to 
obtain  conventional  noise  statistics.  This  feature  permits 
the  comparison  of  selected  3-axis  views  of  noise  with 
standard  statistical  descriptors  of  noise. 
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APPENDIX  B 

FILTERED  IMPULSE  PROCESS 

The  complex  envelope  of  a  noise  process  can  be  described 
by  the  Following  equation 

N(t>  -  j-e, 

n(t)  =  g<t)  +  Z  a,  e  h(t-t,)  <B-1> 

i  =  l 

where  g(t)  is  a  quadrature  Gaussian  process,  h(t)  is  a 
linear  Filter,  a,  is  an  independent,  identically  distributed 
random  amplitude  and  8,  is  a  random  phase  uniFormly 
distributed  over  0  to  2n  .  The  joint  characteristic  Function 
oF  the  inphase  and  quadrature  components  oF  the  above 
process  is 

j(u>,  nc  <t)  +  tit*  n*  (t)  ). 

0(u),  ,o),  )  =  ECe  1  .  (B-2) 

IF  the  Gaussian  noise  is  assumed  independent  oF  the 
impulsive  noise  the  character istic  Function  oF  the  Gaussian 
term  alone  can  be  written 

j(uii  gc  (t)  +  o»»  g*  (t)  ) 

0«  (o),  ,  )  =  ECe  1.  (B-3) 

The  inphase  and  quadrature  components  oF  bandpass  white 
Gaussian  noise  are  independent  so 

2  2 

*No  B(ui  +■  id,  ) /2 

0*  (u),  ,  id,  )  =  e  .  <B-4) 

The  joint  characterist ic  Function  oF  the  ith  pulse  oF  the 


inphase  and  quadrature  impulsive  noise  terms  is 


0i  (<jt  ,o)»  )  =  ECe 


ja(  (u,  (hc  (t-t,  )cosSi  -  h*  (t-t,  )sjn0,  ) 


c 1)1  (he  (t-t,  )sin04  +■  ht  (t-t,  )cos0,  )) 


(B-5) 


Since  the  noise  process  is  narrowband,  its  joint  density 


function  is  circularly  symmetric  and  the  following 
transformat  ion  can  be  applied 


2  2  1/2 
Up  *  (ld|  +  (Hi  ) 


a  =  tan  (u>,  /u.  ) 


(d,  =  <4,.  coso 


Ua  =  (dp  si  no  , 


(B-6) 


This  type  of  transformation  of  a  2-D  Fourier  transform  to  a 
one  dimensional  transform  of  the  envelope  is  known  as  a 
Hankel  Transform  and  the  inversion  formula  is 


p  ( r )  =  r(d,  0  (td,.  )  Jo  (rtdp  >  dur 


( B-7 ) 


where  Jo  is  the  ordinary  Bessel  function  of  the  first  kind, 
order  0.  To  simplify  the  derivation  it  is  assumed  that 
ha (t)  is  zero  (see  comment  before  Eqn.  2-21) .  Now 
transforming  Eqn.  B-5  and  taking  the  expected  value  of  0, 
and  a. 


m  2n 

P  P  j air  a i  hc  (t-ti  )cos(0,  -a) 

0i((dr)=EC  p(0,)p(a, )e  d0 ,  da,  1 


(B-8) 


0  0 


Since  0,  is  uniformly  distributed  over  0  to  2n , 


.v.-O.v-, 

*  a  * 


W 


0  1  (u,  ) 


EC 


p  ( at  )  Jo  ( o)r  a,  hc  (t)  )da,  3. 


<s 


0  ( B  -  9 ) 

If  N(T)  in  Eqn.  B-l  is  a  Poisson  process  with  rate  X  and 
observation  interval  T,  the  t,  's  are  independent  and  uniform 
over  (0,  T)  allowing  the  characteristic  -function  of  the 
impulsive  term  to  be  written  as 

oo  k 

0,  (ID,  )  =  E  0  (u,)  Pr[N  =  k].  ( B - 1 0 ) 

k=0 


Now  using  the  probability  law  for  k  events  of  a  Poisson 
process 


<d  k  -XT  k 

0,  <wr  )  =  E  0  C(ur  )  e  (XT)  /k!  (B-U) 

k=0 


which  sums  to 


p  ( a) 


[Jo  (hi,  ahe  (t))  -  lld-rda 


(B- 12) 


0  0 

0i  (id,  )  =  e 

Using  the  convolution  property  of  zero  order  Hankel 
transforms  CRef.  54]  and  Eqn.  B-7  the  density  function  of 
the  envelope  for  Gaussian  plus  impulsive  noise  is 


<D 


p  C  r  ) 


riu,  Jo  (rill,  )0a  (  0Jr  )0t  (Ur  )du)r  . 


■“00 


(  B  -  1  3 ) 
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APPENDIX  C 

SIGNAL  PLUS  NOISE  DENSITY 

For  the  truncated  Hall  model  with  0  =  3,  the  joint 

density  Function  oF  the  in-phase  and  quadrature  components 
is  given  by: 


2  2 

D  2y  2  2  2 

p(x,y)  =  -  0  <  x  +  y  <  T, 

2  2  2  2  2 


where 


2n(D  - 1 ) ( x  +  y  +  y  ) 
=  0 


2  2  2 
x  +  y  >  Tp 


2  2  1/2 


(C-  1  ) 


D  =  (1  +•  T,  /y  ) 

The  sum  oF  a  random  phase  sinusoidal  signal  oF  amplitude  A 
and  noise  can  be  expressed  in  terms  oF  its  in-phase  and 
quadrature  components  w  and  z  where: 

w  =  AcosCty)  x 

z  =  Asi  n<\|/>  +  y.  (C-2) 


Eqn.  C-2  is  substituted  into  Eqn.  C-l  and  transFormed  to 
polar  coordinates  with 

w  =  rcos(0) 

z  =  rs in(0)  ( C-3 ) 

where  r  represents  the  envelope  and  0  the  phase  oF  the 
signal  plus  noise.  The  joint  probability  density  Function 


p<r,0)  conditioned  on  »y  is 


p(r,0|iy)  =  .  —  1  1  ■  "  . . — 1 

2  2  2  2  2 
2n(D  -l)(r  -  2Arcos(n/  -  0)  +  A  +  y  ) 

2  2  2 
0  <  r  -  2Arcos(\y  -  0)  •+■  A  <  T* 

=  0  elsewhere.  (C-4) 

This  expression  is  valid  -for  all  values  oF  r  between  0  and 
T ,  -  A.  Since  our  principal  interest  is  the  small  signal 

case,  A  <<  T,  ,  this  will  not  aFFect  subsequent  analysis.  By 
integrating  Eqn.  C-4  with  respect  to  0,  the  conditional  mar¬ 
ginal  density  Function  oF  the  envelope,  p(r|ty),  is  obtained; 

n 


p  C  r  |  »y) 


p  C  r  ,0  |  \y)  d0 


0 


<C-5> 


The  integral  is 
nowhere  else  in 
0  and  using  Eqn. 

p  ( r  |  vy  )  =  - 

4 


a  Function  only  oF  cosCiy  -  0),  iy  appears 

the  integral  .  ThereFore  by  setting  ct  =  >y  - 

2.554-3  in  CReF.  54D  gives 

2  2  2  2  2 
Cr  +  A  •+>  y  )  2y  rD 

22  22  22442 


(r  -  2A  r  +  2A  y  *  2y  r  +  A  +  y  ) (D  -  1) 


(C-6) 


Since  ly  is  assumed  independent  From  sample  to  sample  and 
uniFormly  distributed  over  2n  , 

2n 


p  ( r ) 


p(r|\y)p(»y)  d\y 


(C-7) 


p(r) 


2  2  2  2  2 
( r  +  A  +  y  J  2y  rD 


4  22  22  22  442 

( r  -  2A  r  •+•  2A  y  +■  2y  r  2-*-  A  +  y  )  ( D  -  1) 


By  a  similar  argument,  the  density  -function  -for  a 
phase  sinusoidal  signal  in  narrowband  noise  with  a 
distribution  is 

4  2  2 

2y  r ( b  +  2a  ) 
p(r)  =  1  1  - 

2  2  5/2 

(a  -  b  ) 

where 

2  2  2 
a  =  r  +■  A  +  y 

and 

b  =  -2Ar. 
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APPENDIX  D 


SPECTRUM  OF  A  RENEWAL  PROCESS 

The  complex  envelope  representation  oF  Filtered  impulse 
noise  on  a  Finite  interval  T  is  given  by 


N  ( T )  -j0, 

n(t)  -  E  a,  e  h(t-t!  )  ; 

1  =  1 


0  <  t  <  T. 
(D-  1  ) 

N(T)  is  a  random  unit  counting  process  denoting  the  number 
oF  impulses  in  the  interval,  the  -Ctt  >  are  the  random  arrival 


envelope  impulse  response  h<t) . 

The  envelope  squared  oF  n(t>  is 

2  N(T)  2 

E  <t>  =  n(t)n"  (t)  =  E  a,  h ( t-t,  )  h«  ( t-t,  ) 

i  =  l 

( D-2 ) 

where  we  have  assumed  terms  oF  the  Form  h(t-t, ) h*  (t-tk )  For 
diFFerent  i  and  k  are  negligible.  We  estimate  the  power 
spectral  density  (PSD)  oF  the  envelope  by  averaging 
magnitude  squared,  length  T,  Fourier  transForms  oF  E* (t). 
The  mean  value  oF  the  estimate  is 


S  (u>,  T)  =  ECE,  <u»>£.  •  (w>  3/T 


(D-3) 
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ii 


where 


£a  C(0,T) 


T 

2  -jut  N ( T )  2  -jut, 

E  Ct)e  dt  =  £  a(  e  Ha  Cu) 


and 


H.  (u>) 


-  jut 

h(t)h« (t)e  dt 


CD 

is  •the  Fourier  transform  of  the  magnitude  squared  impulse 
response.  We  have  assumed  the  impulse  response  is  much 
shorter  than  the  observation  interval  T. 

The  mean  of  the  estimate  becomes 


S  ( o> ,  T ) 


NCT)  N  ( T )  2  2  -j(i>Ct,-tk>  2 

EC  £  £  a»  ak  e  1  |  H,  (id)  |  /T. 

U= 1  1=1 

( D-4  ) 


This  is  a  compact  and  general  expression  For  the  PSD  oF  a 
truncated  Fi ltered  point  process.  It  is  in  terms  oF  a 
summation  oF  the  characterist ic  Functions  oF  the 
interarrival  time  distributions  between  all  combinations  oF 
pairs  oF  points.  (The  term  incorporat ing  the  Frequency 
response  oF  the  Filter  will  be  set  to  one  For  the  rest  oF 
the  analysis.)  Conditioning  on  NCT) 

<d  N  N  2  2  -ju(t,  -tk  ) 

S  ( to ,  T )  =  £  EC  E  E  a,  a*  e  IPrCN  =  1J/T. 

1=2  k= 1  i=l 

(D-5) 

Plotting  the  diFFerence  k-j  indicates  a  diFFerent  summation 


over  m  and  n  shown  below 


2  -1 
1  0 

0  1 


Let  m  =  i  -  k  and  n  =  i  +  k  such  that  t.  =  t,  -  t„  is  the 
mth  interarrival  time.  Then  Eqn.  D-5  can  be  written  as 


_  ®  N-l-n  2  -jut. 

S ( u , T)  =  E  E-C  E  E  Ag  e  +  NA« 

1=2  n= 1  m=l 

N  - 1  n  2  .  jut. 

+  E  E  A.  e  >PrCN  =  1 1/T. 

n=l  m=l 


(D-6) 


On  -Further  rearrangement 

_  _  <n  2  N- 1  n  jut. 

S(u,T)=  NA, /T  +  E  A,  E-C  E  E  e 

1=2  n=l  m=l 


-  jut. 


>PrCN  =  11/T 


(D-7) 


and  conditioning  the  interarrival  time  distribution  over  the 
number  o-P  pulses  N  in  the  interval  T 

2  «  1-1  n 

S(u,T)  =  NA« /T  +  A,  E  E  E  C0. Cju|N=l)  + 

1=2  n= 1  m=l 

0. ■  C ju| N=1 ) IPrCN  =  11/T 


where 


jut. 

0. ( ju| N= 1 )  =  ECe  1 


is  the  characteristic  Function  oF  the  m* h  interpulse  spacing 


given  exactly  1  pulses  in  the  interval 


assumption  is  now  made;  that  the  mth  interpulse  spacing  is 


The  Fo 1 1 ow l ng 
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■the  sum  of  m  statistical  ly  independent  interpulse  sep¬ 
arations,  and  is  independent  of  1,  the  number  of  pulses  in 
the  interval.  This  assumption  will  give  an  asymptotically 
correct  result  -For  the  PSD  as  the  observation  interval  T, 
becomes  much  greater  than  p,  the  mean  value  of  the 
interarrival  times. 

With  the  above  assumption 


0.  ( ju)|  N=1  )  =0  (ju) 


(D-9) 


where  0(ju»)  is  the  character i st i c  -function  of  the  inter¬ 
pulse  separations.  Eqn.  D-8  then  sums  to 


_  _  2  ®  1-1  0  ( j  u) )  -0  (  j  u) ) 

S  ( <ii ,  T  )  =  NA«  /T  +  Aa  E  EC  -  + 

1=2  n=  1  l-0(joj) 


0*  ( j(D)  -0*  (jin) 


1  -0"  ( ju>) 


UPrCIM  *  13/T 


CD- 10) 


and  summing  again 


_  _  2®  10  (  j  (D )  0  C  ju»)  0  (  j  u» ) 

S  ( (u ,  T )  =  NA«  /T  +  A,  E  -C -  -  -  +  - 

1=2  2  2 
<l-0(j<u>)  (l-0(jti>))  (1-0  Cju)  ) 


10-  (;in) 


0-  ( jU)  ) 


0-  Cju) 


2  2 

(  1-0*  <  jW>  >  (l-0-(jU))  (  1  —0  -  (  j  (i)  )  ) 


>PrCN  =  11/T 


(D-  1  1  ) 


Rearrang i ng 


_  _  2  _  ®  0  ( ju)  0*  (  j  a) ) 

S(w,T)  =  NA«  /T  +  A, /T<N  -  PrCN=l!)>  E  < -  +  - 

1=2 

(  1  -0(ju)  )  (  1  -0  (  jui)  ) 


2  0Cj(j) 

-  A, /T  -Cl  -  PrCN=03-PrCN=l>  < - 

2 

(  1  -0(jiii)  ) 


2+1  1+1 

2  cd  0(joj)  0*  (ju) 

+  A, /T  E  Pr  CN  =  13  < -  +  -  >. 

1=2  2  2 

(l-0(jw))  Cl-0"  <j(D))  (D-12) 

Since  PrCN  =  13  is  less  than  1  and  |0(<u)|  <  1  CReF. 

41 :p.  1153  for  any  distribution,  then  the  absolute 

convergence  oF  the  third  term  is  guaranteed.  For  an 

equilibrium  renewal  process 

N  =  ECN (t ) 3  =  T/p  (D-12) 

where  p  is  the  mean  time  between  renewals  CReF.  32 :  p .  46  3. 

Now  letting  T,  the  observation  interval ,  go  to  infinity  the 

asymptotic  PSD  is 


_ 

_ 

Aa 

A* 

S(«d) 

=  lim  S(ui,T)  = 
T=>« 

M 

C  -  +  2 

2 

Aa 

ReC 

This  is 

the  expression 

For 

the  PSD 

oF 

i ncrement 

point  process 

i  n 

terms  oF 

the 

Function  oP  the  interarrival  distribution, 
result  was  developed  by  Cox  and  Miller  by 
renewal  intensity  Function  CReF.  553. 


0  <  j(ll) 

-  3  >  . 

1-0C  j(i») 

CD-13) 
an  independent 
character istic 
A  similar  result 
considering  the 
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APPENDIX  E 


SPECTRUM  OF  A  BRANCHING  RENEWAL  POINT  PROCESS 

The  complex  envelope  representation  of  Filtered  impulse 
noise  arising  From  a  branching  renewal  process  on  a  Finite 
interval  T  is  given  by 

M(T)  N.  -j0.„ 

n<t)  =  E  E  a.  „  e  h(t-t.„-T,  )  (E-l) 

m=l  n=l 

where  T.  is  the  beginning  oF  the  mth  main  interval,  N.  is 

the  number  oF  impulses  in  the  mth  interval ,  arriving  at 

times  -Ct.  „  >  aFter  its  onset  and  M(T)  is  the  unit  counting 

process  deFining  the  number  oF  main  events  in  the 

observation  interval  T.  The  pulse  amplitudes  <a. „>  and 

phases  <0.„>  are  statistically  independent  random  variables 

and  the  {0, „ >  are  taken  to  be  uniForm  on  <0,2n>.  The  net 

Filtering  eFFects  on  the  impulses  are  accounted  For  by  the 

time  invariant  complex  envelope  impulse  response,  h(t). 

The  envelope  squared  oF  the  process  is 

2  M<T)  N.  2  2 

E  (t)  =  n(t)n" (t)  =  E  E  a  . „ | h (t-t„ „ + T. ) |  (E-2) 

m= 1  n=l 

and  using  the  same  assumptions  as  Appendix  D  the  mean  value 
oF  the  Bartlett  estimate  oF  the  PSD  is 

2 

S  ( (u  ,  T )  =  EC  I  E,  (w,  T)  I  1/T  ( E-3 ) 

where  Et  (w,T)  is  the  length  T  Fourier  Trar.sForm  oF  Ea  Ct)  . 
Spec i Fi ca 1 1 y , 
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j  »  »  *  -  _  ^ _  r H  w^jv. 


-  .  •  ,  •r-  r 


M(T)  N.  2  -jut. „  -juT. 

£»  <(d,T)=  E  E  a  . .  e  e  Ha  (u) 

m= 1  n=l 


(E-4) 


with  H*  (u) 


T 

2  -jut 

| h<t) |  e  dt ,  the  Fourier  transform  of  the 


magnitude  squared  impulse  response.  To  find  the  mean  of  the 

estimate  we  must  evaluate 

2  M(T)  N.  M(T)  N, 

SCu,T)  =  EC  | Ha  (u)  |  E  E  E  E 

m=l  n=l  1=1  k=l 

2  2  —  j  u  (  t.  .  —  1 1  k  )  -  JU  (  T.  — Ti  ) 

a.  „  ai  .  e  e  1/T 

(E-5) 


In  order  to  simplify  this  expression,  we  make  use  of 
assumptions  1,  2  and  4  from  section  3D;  (1)  EC  a*  .  „  aa  ,  „  3  = 
a*  ,  which  will  be  normalized  to  one  for  simplicity,  (2)  N. 
and  N,  equal  either  N1  or  N2  alternating  between  the  two 
values  and  <3)  f(t. „)  and  f (tj „ ) ,  the  interarrival  time 
distributions,  alternate  between  the  negative  and  positive 
phase  distribution  parameters.  Furthermore,  let  M  be  an 
even  number  of  main  process  points  occurring  at  intervals, 
T./2,  of  a  fundamental  frequency  (in  our  case  f.  =  u. /2n  = 
1/T.  =  60  Hz,  the  power-line  frequency)  and  let  M2  =  M/2. 
This  assumption  will  greatly  simplify  the  calculations  and 
for  M  greater  than  10,  not  significantly  reduce  the  accuracy 


of  the  result. 


Now  looking 


only  at  the  factors  involved 


1  n  the 


summation  in  Eqn.  E-5,  six  separate  terms  can  be  identified: 


Term  1  -  The  inter-pulse  group  summations  between  pulses  in 

pulse  groups  with  N1  pulses 

N 1  N 1  ju>(t„  -t„  ) 

=  M2  E  E  ECe  1 

n= 1  k=l 

(E-6a) 

Term  2  -  The  inter-pulse  group  summations  between  pulses  in 

pulse  groups  with  N2  pulses 

N2  N2  jai(t„  -tk  ) 

=  M2  E  E  ECe  1 

n= 1  k= 1 

(E-Sb  "i 


Term 

3 

-  The 

i  ntra-pu 1 se 

group 

summat i ons 

between  pulses 

i  n 

groups 

with  N1 

pu 1 ses 

M2  M2 

-ju(m-l)T,, 

N1 

N 1 

jdlto 

-  ju)tk 

= 

E  E 

e 

E 

E 

ECe 

lECe  1 

m= 1  1=1 

n=  1 

k=  1 

m<  >  1 

CE- 

6c ) 

Term 

4 

-  The  intra-pulse  group  summations  between  pulses 

i  n 

groups 

with  N2 

pulses 

M2  M2 

-  jo)  Cm- 1 )  T0 

N2 

N2 

j  tiitn 

-  j  wtk 

= 

<  E  E 

e 

E 

E 

ECe 

]ECe  ] 

m=l  1=1 

n=  1 

k=  1 

m<  >  1 

<E- 

6d ) 

T  erm 

5 

-  The 

l ntra-pu 1 se 

group 

summat i ons 

between  pulses 

l  n 

groups 

with  N1  pulses  per 

group 

to 

pu 1 ses 

in  groups  with 

N2 

Terms  5  and  6  are  complex  conjugates  and  can  be  combined  to 


form  one  term 

cos(M2  2nw/w0  )  -  1 

Term  5a  =  2cos  (w  ( T,  /2-T0  )  )  <  —  - . -  >4/(0,  ,N1,0*  ,N2) 

cos  (2na »/<u.  )  -  1 

(E-7e> 

where  0,  and  0a  are  the  character i st i c  -functions  associated 
with  the  interarrival  times  for  each  phase  and  N1  and  N2  are 
the  number  of  renewals  associated  with  each  phase.  The 

functions  5,  fi,  and  14/  are  defined  below; 


2  N 

N  -  n;0(j(d)  I  0(jo»)  (1  -  0i.ju))  ) 

§(0,N)  =  -  -  2  Re-C -  >, 


I  1  -  0<jw)  | 


fK0,N) 


(1  -  0 C jw> ) 


N+l  2 

|0(j«d)  -  0(jw)  | 


(E-8) 


I  1  -  0(ju»)  | 


CE-9) 


and 


14/(0,  ,  N 1  ,0a  ,N2)  = 


Thu« 


Nl+1  N2+1 

(0,  <J(J)  -  0,  <JU>)  )  (0.  (-JU>)  -  0a  ("JO))  > 


<1-0,  <  j(D>  )  (1  -  0,  <-j(D)  ) 


S  ( o) ,  T )  =  -CECH,  <u)3/T>  L  Term(n). 


(E- 10) 
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APPENDIX  F 


NON-HOMOGENEOUS  POISSON  PROCESS 


The  Following  derivation  is  for  -the  mean  and  variance  of 
a  compound  periodic  non-homogeneous  filtered  Poisson  point 
process.  This  proof  is  a  generalization  of  a  proof 
presented  in  Papoulis  CRef.  41:pp.  382-3831.  The  expression 
for  such  a  process  on  an  interval  0  to  T  is  given  by 


NCT) 

n(t)  =  E  a, h(t-ti ) , 
i  =  l 


(F-  1  ) 


where  N(T)  is  a  unit  counting  process  with  periodic  rate 
parameter  ACt) ,  a,  is  the  weighting  of  the  ith  point  and  is 
independent  and  identically  distributed  (I  ID)  from  point  to 
point  and  h(t)  is  the  time  invariant  impulse  response  of  the 
f i Iter . 

The  time  axis  is  divided  into  consecutive  intervals  Ij 
of  length  aoi  where  aoc  =  otj  -  Oj  .  i  .  The  number  of  jumps  in 
the  counting  process  in  the  j* h  interval  I,  is  given  by  aih,  . 
If  Aa  is  sufficiently  small,  then  the  contribution  to  the 
total  noise  process  from  the  j* *  interval  is 


Anj  *  E  a(  h Ct-Oj  )  , 
i  =  l 


( F  -  2  ) 


where  arrij  is  a  Poisson  random  variable  with  its  rate 


parameter  approx imately  equal  to 


A  ( a  ,  )  a  oc . 


(F-3) 


u 


The  moment  -function  of  the  j*h  interval  is  then 

Aftlj 

s  E  a,  h ( t-o j  ) 
i  =  l 

a0j  (s)  »  ECe  3.  CF-4) 

Conditioning  on  adIj  , 

Alltj 

s  E  a,  h ( t-cxj  ) 

1  =  1 

a0j  (s)  *  ECe  | Afflj  =  k3  PrCAmj  =  k3  CF-5) 

and 

k  sa,  h ( t-Oj  ) 

A0J  (s)  *  EC  n  e  3  PrCAHtj  =  k3.  (F-6) 

1  =  1 

Since  the  a,  's  are  independent  -from  pulse  to  pulse 

k  sah(t-otj  ) 

as j  <s)  *  n  ECe  3  PrCAmj  =  k3,  CF-7) 

1  =  1 

where  the  expectation  is  with  respect  to  the  random 
amplitude  a.  Using  the  Poisson  probability  law  and  recog¬ 
nizing  the  series  expansion  For  the  exponential  gives 


sah (t-cxj  ) 

XCajiAOcCECe  3-13- 

a0j  (s)  »  e  .  (F-8) 

Since 

CO 

n(t)  =  E  An, (t),  (F-9) 

i  =  -m 

a  sum  oF  independent  random  variables,  then 

00 

0„  Cs)  =  n  A0 ,  (S)  CF-  10) 

i  =  -® 

using  the  convolution  property  oF  the  moment  generating 
Function.  The  cumulant  generating  Function  is  deFined  as 


APPENDIX  G 


AUTOCORRELATION  OF  CORONA  NOISE 


The  mean  value  oF  the  cyclostationary  process  defined  by 
assumptions  1  and  2  in  section  4D  is  given  by  Eqn.  4-12. 
Substituting  Eqn.  4-20  gives 


ECeCt) 3 


ECal  J 
0 


2 

Nt  1/2  -as 

E  a„  cos  (u«  n (t-s )  )  Ca /it  )  e  ds.  (G-l) 
n=0 


Expanding  the  cosine  and  integrating 

2  2 
-(Do  n 


NT  4a 

ECe<t)3  =  (fln/4)  E  C  e  a„cos((D0nt)  + 

n=0 


2  2 

(Do  n  oo  1  -(Do  n  k - 1 

+  _____  E  -  (  ■  )  a„  sin (u)o  nt)  3 

1/2  k=l  ( 2k- 1 >  »  !  2a 

2<an  ) 

(G-2) 

1/2 

where  <fln/2)  =  ECa3  and  (2k-l)M  =  1  3  5...2k-l  CReF.  54: 
Eqn.  3.8973.  The  sine  term  is  in  quadrature  with  the 
dominant  cosine  term  of  Eqn.  G-2  and  represents  a  phase 
shift  in  the  mean  oF  the  process  relative  to  the  rate 
Function.  This  is  due  to  the  Filtering  eFFect  where  the 
impulse  response  oF  the  Filter  persists,  thus  causing  a  lag 
in  the  mean  value  oF  the  process.  For  the  values  oF  a  and 


(Do  considered  in  this  work,  the  sine  term  in  Eqn.  G-2  will 


be  negligible  relative  to  the  cosine  term.  Physically,  this 


means  that  -for  filter  bandwidths  much  greater  than  the 
fundamental  frequency  of  the  corona  noise,  the  phase  shift 
will  be  insignificant. 

Similarly  the  variance  of  the  envelope  when  considered 
as  a  cyclostationary  process  is 

2  2 
-u>o  n 


VarCe<t):i  *  2fl 


NT *  1  1/2  8a 

E  C  -  Ca/2n>  e 

n=0  2 


a„  cos  ( (d0  nt )  3  . 

(G-3) 


Using  Eqn.  4-14  the  covariance  of  the  assumed  process  is 


Cov  Ct,  ,  to  )  =  2fl 


2  2 

P  Nt  -a(t,  -s)  -aCt»  -s) 

E  a„  cos((do  ns)  Ca/n  )e  e  ds 

n=0 

max (t, , t* )  ( G-4 ) 


where  after  rearranging  and  completing  the  square 


-a  <t,  -t.  )  /2 

Cov(t, ,tj )  =  (2fl(a/n)e  )  x 


CO  2 

'  Nt  -2a  (s-  Ct,  +ta  )  /2 > 

E  a„  cos  ( (do  ns )  e  ds. 

n=0 

maxCt, ,t* )  (G-5) 


Lett i ng 

u  =  (s  -  Ct,  +  t,  )/2),  CG-6) 

the  covariance  is 


Cov  (t,  ,-t*  ) 


-a  ( 1 1  -ta  ) 

(  (a/n )  213  e 


x 


CD  2 

'  Nt  -2au 

£  a„  cos  ( u)o  n  ( u  +  <t,  +ta  ) /2 )  e  du 
n=0 

|t,  -t.  |  (G-7) 

which  cannot  be  evaluated  in  closed  form. 

The  stationary  autocorre  1  at i on  -Function  is  Found  by 
using  Eqns.  G-2  and  G-7  and  the  deFinitions  given  by  Eqns. 
4-15  and  4-19.  The  term  due  to  the  covariance  in  the  auto¬ 
correlation  Function  is 


2 

-ct'j  /2 

2l3(ot/n)e  /T0 


Nt 

E  Ca„  cos  <u)0  n  (2t+T  )  ) 
n=0 


-2au 

cosCuo  nu)e  du 


0 


1^1 


<s  2 

I"  -2au 

sintuio  nu)e  du3dt. 

| t |  (G-8) 

The  sinusoidal  terms  will  integrate  to  zero  with  the 

exception  oF  the  a0  cosine  term  and  the  expression  simp- 

1 l Fi es  to 


-oct  /2 

213  ( a  /  n  )  e  a0 


-2au 

e  du . 


M 


( G-9 ) 

This  term  will  be  dominated  by  the  First  exponential  and  is 


approximated  by 


APPENDIX  H 


LOCALLY  OPTIMUM  ENERGY  RECEIVER 


The  observations  consist  o-P  a  sequence  o-P  N  complex 


samples.  The  hypo-theses  are: 


H„  :  Xi  =  ni 


versus 


H,  :  x,  =  nt  +  A 


i  =  1  ,  .  .  .  ,  N 


(H-l ) 


where  nt  are  an  independent,  identically  distributed  Cl  ID) 
complex  noise  samples  with  a  joint  inphase  and  quadrature 
density  p(ne ,n,  ).  The  sequence  s,  is  a  complex,  zero  mean 
signal  sequence  with  a  known  variance,  <x,  ,  and  A  is  a  real, 
positive  number.  The  1/2  power  o-P  A  is  chosen  to  facilitate 
the  derivation. 

The  generalized  likelihood  ratio  CRe-P.  6:p.  585,213  to 
test  the  hypotheses  described  by  Eqn.  H-l  is 


A  <  x  ,  A )  =  n 
1  =  1 


1/2  1/2 
pCxc  i  -A  5Ci  ,x,  ,  -A  s«  ,  ) 

p(xc  i  ,  x*  ,  ) 


p  C  Sc  i  ,  s*  i  )  dsc  i  ds* 


( H-2) 


Applying  Eqn.  5-8,  and  d  i  f-Perent  i  at  i  ng  with  respect  to  A, 


tL  o  (x  ,  A)  =  L 
i  =  1 


1/2  1/2 

-SciOpCxei  -A  Sci  ,X,,-A  s«i>/2 


uxc  I  p(xc  i  , X*  ,  )  A 


-S«  !  Op<Xe  »  -A  5(1  ,X|  ,  -A  Sat>/2 


1/2 

■ax*  ,  p(xe  i  ,  X,  ,  )  A 


3  p<Sc 


1S1  i  )  dSe  t  ds*  i 


A=>0. 


(H-3) 


Letting  A  go  to  0  and  simultaneously  taking  the  expectation 
oP  se  i  and  s*  ,  requires  L'hospital's  rule  to  evaluate  the 
expression  EReP.  563.  APter  taking  derivatives  the  test 
statistic  is 


tlo  <  x) 


<D  CD 


N 

=  £  C 

1  =  1 

-CD  -CD 


2  2  1/2  1/2 
Sc  i  10  p(xc  i  ~A  Sc  i  i  x»  i  -A  s*  i  )  /2 

2 

axc  i  p  (  xc  i  ,  x*  i  ) 


2  1/2  1/2 
Sc  i  s*  i  a  p(xe  i  -A  Sc  t  ,  x*  t  -  A  s*  (  ) 


axc  i  TJX*  i  p(xc  t  ,x*  i  ) 


2  2  1/2  1/2 
s*  i  a  p  Cxc  i  -A  sc  i  ,  x«  i  -  A  s*  ,  )  /2 

“ “ — — ~ — — — — — — — — •  ]  p  ( Sc  i  ,  s*  i  )  dsc  i  ds* 

2 

ax,  ,  p  (xc  i  f  x,  ,  ) 


A =>0. 
(  H-4  ) 


In  general,  the  inphase  and  quadrature  noise  components  will 
be  circularly  symmetric  even  though  not  necessarily 
statistically  independent  [ReP.  41 :p.  1333  such  that 

2  2  1/2 


with 


P(r )  =p((xe  +  x*  ) 

2  2  1/2 


CH-5) 


r  =  (xe  +  x. 


( H-6 ) 


LIST  OF  REFERENCES 


Cummins,  E.J.,  Jauregui ,  S.,  and  Vincent,  W.R.,  "Time- 
and  Frequency-Domain  Characteristics  oF  Man-Made  Radio 
Noise  AFFecting  HF-Commun i cat i ons  Sites,"  IEEE  Trans. 
on  EMC,  v.  EMC-21,  pp.  182-189,  August  1979. 

Herman,  J.R.,  "Survey  oF  Man-Made  Radio  Noise," 

Progress _ in  Radio  Science.  International  Union  oF  Radio 

Scientists,  Brussels,  Belgium,  pp .  315-408,  1970. 

StanFord  Electronics  Laboratories  Technical  Report  No. 

7 05 0 - 7 H ,  A  New  Model  For  Impulsive  Phenomena; 
Application  to  Atmospheric-Noise  Communicati on 
Channe 1 s .  by  H.M.  Hall,  August  1966. 

De  ReFFye ,  J.,  General  Modelization  oF  Impulsive  Noise, 
Inter.  ConF.  on  Commun. ,  1982,  v.  1.  pp.  2H. 3/1-5,  June 

1982. 

Modestino,  J.W.,  Jung,  K.Y.  and  Matis,  K.R.,  Modeling, 
Analysis  and  Simulation  oF  Receiver  PerFormance  l n 
Impulsive  Noise.  GLOBECOM  1983,  v.  3,  pp .  1598-1605, 

December  1983. 

Urkowitz ,  H.,  Signal  Theory  and  Random  Processe s , 

Artech  House,  1983. 

Brigham,  E.O.,  The  Fast  Fourier  TransForm.  Prentice- 
Hall,  Inc.,  p.  230,  1974. 

Kailath,  T.,  "The  Complex  Envelope  oF  White  Noise," 

IEEE  Trans,  on  InF.  Theory,  v.  IT- 12,  pp .  397-398,  July 
1966. 

Di Franco,  J.V.  and  Rubin,  W.L.,  Radar  Detection. 
Prentice -Hal  1 ,  Inc.,  1968. 

Beckmann,  P,  "Amplitude  Distribution  oF  Atmospheric 
Radio  Noise,"  Radio  Science,  y.  68D ,  p.  726,  June  1964. 

Electric  Power  Research  Institute  Report  EL/EM  -4290, 

v .  1 ,  Harmonics  and  Electrical  Noise _ in  Distribution 

Systems .  by  W.R.  Vincent,  April  1985. 

Myers,  H.A.,  "Industrial  Equipment  Spectrum 
Signatures,"  IEEE  Trans,  on  RFI.  y.  RFI-5,  pp.  30-42, 
March  1963. 


Crichlow,  W.  Q.  ,  and  others,  "Determination  of  the 
Ampl  ltude-Probabi  1  ity  Distribution  o-f  Atmospheric  Radio 
Noise  From  Statistical  Moments,"  J.  o-f  Res,  o-f  NBS,  v. 
61D,  January-February  1960. 

Likhter,  J.  I.,  "Some  Statistical  Properties  of 
Atmospherics,"  Radiotekh  i  Electron,  v.  1,  pp .  1295- 

1302,  1956. 

Ibukan,  0.,  "Structural  Aspects  of  Atmospheric  Radio 
Noise  in  the  Tropics,"  Proc.  IEEE,  v.  54,  pp .  361-367, 
March  1966. 

Spaulding,  A.  D.,  Roubique,  C.  J.,  and  Crichlow,  W.Q., 
"Conversion  of  the  Amplitude-Probability  Distribution 
Function  for  Atmospheric  Radio  Noise  from  One  Bandwidth 

to  Another,"  J.  of  Res. _ of  NBS.  v.  66D ,  pp .  713-720, 

November-December  1962. 

Horner,  F.  and  Harwood,  J.,  "Investigation  of 
Atmospheric  Radio  Noise  at  Very  Low  Frequencies,"  Proc. 
IRE,  v.  103,  pt  B,  pp.  743-751,  November  1956. 

Schonhoff,  T.A.,  Giordano,  A. A.,  and  Huntoon ,  Z., 
Analytical  Representation  of  Atmospheric  Noise 
Distribution  Constrained  in  V„ .  Inter.  Conf.  on 
Commun.,  Chicago,  Illinois,  June  1977. 

Katz,  A.H.,  Meltz,  G.,  Mullen,  J.A.,  "Probability 
Density  Distribution  of  Wideband  2-6  MHz  Terrestrial 
Rad l o  No i se , "  IEEE  Trans,  on  Antennas  and  Propaq . ,  v . 
AP-27,  pp .  107-111,  January  1979. 

Lord,  R.D.,  "The  use  of  the  Hankel  Transform  in 
Statistics,"  Bipmetr ika ,  v.  41,  pp.  44-55,  1954. 

Modestino,  J.W.,  and  Ningo,  A.Y.,  "Detection  of  Weak 
Signals  in  Narrowband  Non-Gaussian  Noise,"  IEEE  Trans. 
on  Inf.  Theory,  v.  IT-25,  pp .  592-600,  September  1979. 

Furutsu ,  K.  and  Ishida,  T.  "On  the  Theory  of  Amplitude 

Distribution  of  Impulsive  Random  Noise,"  J. _ Appl .  Phy s. 

( Japan ) ,  v.  23,  pp .  1206-1221,  July  1960. 

Giordano,  A. A.  and  Haber,  F.,  "Modeling  of  Atmospheric 
Noise",  Radio  Science,  v.  7,  pp.  1011-1023,  November 


Middle-ton,  D.,  "Statistical -Physical  Models  oF 
Electromagnetic  Interference, "  IEEE  Trans,  on  EMC,  v. 
EMC- 1 9 ,  pp.  106-127,  August  1977. 


Geneva,  1964. 


Huntoon,  2.M.  and  Giordano,  A. A.,  Energy  Detection in 
Impulsive  Noise.  Inter.  ConF.  on  Commun.  ,  pp.  07.1/1-4, 
June  1979. 

Herman,  J.R.,  and  others,  "Measurement  and  Statistical 
Analysis  o-P  Wideband  MF  Atmospheric  Radio  Noise,"  Radio 
Science .  v.  21,  pp.  25-46,  January-February  1986. 

Skomal,  E.N.,  Man-Made  Radio  Noise.  Van  Nostrand 
Re inhold  Company,  NY,  1978. 

Pakala,  W.E.  and  Chartier,  V.L.,  "Radio  Wave 
Measurements  on  Overhead  Power  Lines  From  2.5  to  800 
kV,"  IEEE  Trans,  on  Power  Apparatus  and  Systems,  v. 
PAS-90,  p.  1155,  May-June  1971. 

IBM  Technical  Report  RC-2878,  Asymptotic  Properties  oF 
Branching  Renewal  Processes,  by  P.A.W.  Lewis,  pp .  2-5, 
May  1970. 

Moose,  P.M.,  and  0' Dwyer,  J.M.,  "A  Model  For  Impulsive 
Power  Line  Radio  InterFerence  due  to  Gap  Type 
Discharges,"  IEEE  Trans,  on  EMC,  to  appear  November, 
1986. 

Cox,  D.R.,  Renewal  Theory.  Metheun  and  Co.  Ltd.,  1962. 

La r son ,  H . J . ,  Introduction  to  Probability  Theory  and 
Statistical  InFerence.  John  Wiley,  Inc.,  1974. 

Arai,  K. ,  Jani schewsky j ,  W.  and  Miguchi,  N. ,  "Micro-gap 
Discharge  Phenomena  and  Television  InterFerence, "  IEEE 
Trans,  on  Power  Apparatus  and  Systems ,  v .  PAS -104,  No . 
1,  pp.  224-224,  January  1985. 

Ross,  S.M.,  Introduction  to  Probability  Models. 

Academic  Press,  Inc.,  p.  448,  1985. 

Timperley,  J.E.,  "Incipient  Fault  IdentiFication 
through  Neutral  RF  Monitoring  oF  Large  Rotating 
Mach i nes,"  IEEE  Trans,  on  Power  Apparatus  and  Systems , 
v.  PAS  102,  March  1983.  . 


Loeb,  L.B.,  Electrical  Coronas.  University  of 
California  Press,  p.  572,  1965. 

Olsen,  R.G.,  "Radio  Noise  Fields  Generated  by  Corona 
Streamers  on  a  Power  Line,”  Radio  Science,  v.  18,  pp. 
399-408,  May-June  1983. 

Weiss,  G. ,  "Shot  Noise  Models  for  the  Generation  of 
Synthetic  Streamf low  Data,"  Water  Resources  Research, 
v.  13,  pp.  101-108,  February  1977. 

Gardner,  W.A.  and  Franks,  L.E.,  "Character i zat i on  of 
Cyclostationary  Random  Signal  Processes,"  IEEE  Trans, 
on  Inf.  Theory,  v.  IT-21,  pp .  4-14,  January  1975. 

Papou 1 i s ,  A . ,  Probability,  Random  Variables  and 
Stochastic  Processes.  McGraw-Hill,  Inc.,  1985. 

Urkowitz,  H.  "Energy  Detection  of  Unknown  Deterministic 
Signals,"  Proc.  of  IEEE,  v.  55,  April  1967. 

Cook,  C.E.,  and  Bernfeld,  M. ,  Radar  Signals.  An 
Introduction  to  Theory  and  Applications ,  Academic 
Press,  p.  164,  1967. 

Williams,  J.R.  and  Ricker,  G.G.,  "Signal  Detectability 
Performance  of  Optimum  Fourier  Receivers,"  IEEE  Trans. 
on  Audio  and  Electroacoustics,  v.  AU-20,  October  1972. 

Abramowitz,  M.,  and  Stegun,  I.,  Handbook  of 

Mathemat i ca 1  Functions .  US  Government  Printing  Office, 

Washington,  DC,  1964. 

Miller,  J.H.,  and  Thomas,  J.B.,  "Detectors  for 
Discrete-Time  Signals  in  Non-Gaussian  Noise,"  IEEE 
Trans,  on  Inf,  Theory,  v.  IT-18,  pp.  241-250,  March 
1972. 

Rudnick,  P.,  "Likelihood  Detection  of  Small  Signals  in 
Stationary  Noise,"  J .  Add  1 .  Phys . .  v.  32  pp.  140-143, 
February  1961. 

Capon,  J.  "On  the  Asymptotic  Efficiency  of  Local ly- 
Opt i mum  Detectors , "  IRE  Trans,  on  Inf.  Theory ,  v .  I T - 7 , 
pp.  67-71,  April  1961. 

Evans,  J.E.  and  Griffiths,  A.S.,  "Design  of  a  Sanguine 
Noise  Process  Based  Upon  World-Wide  Extremely  Low 
Frequency  (ELF)  Recordings,"  IEEE  Trans,  on  Commun. , 
COM-22,  pp.  528-539,  April  1974. 


v 


-  k  »  fa  m  ^  k  f.. 


M  *  V  ”  k  - 


•  *  k  •  k  * 


v-x  »  ,  -  „  v  r  i  ■»,  »’  j*  x  ’  v  ^  ’ 


50.  Dillard,  Q.M.,  "Mean-Level  Detection  o-F  Non-Fluctuating 
Signals,"  IEEE  Trans,  on  Aero,  and  Elec.  Sys. .  v.  AES- 
10,  pp.  795-799,  November  1974. 

51.  Wozencraft,  J.M.  and  Jacobs,  I.R.,  Principles  o-F 
Commun j cat i on  Eng i neer i ng .  John  Wiley,  Inc.,  p.  459, 
1965. 

52.  SRI  Technical  Memorandum  3564-061583,  Description  oF 

I nstrumentat i on .  by  W.R.  Vincent,  pp.  2-5,  June  1983. 

53.  Papoulis,  A.,  Systems  and  Transforms  with  Applications 
in  Optics,  p.  143,  McGraw-Hill,  Inc.,  1968. 

54.  Gradshteyn,  I.S.  and  Ryzhik,  I.M.,  Tables  oF  Integrals. 

Series  and  Products .  Academic  Press ,  1965. 

55 .  Cox ,  D .  R .  and  Miller,  H.D.,  The  Theory  o-F  Stochastic 
Processes .  John  Wiley,  Inc.,  p.  359,  1965. 

56.  Poor,  H.V.,  and  Thomas,  J.B.  "Locally  Optimum  Detection 

o-F  Discrete-Time  Stochastic  Signals  in  Non-Gaussian 
Noise,"  J.  Acoust.  Soc.  Am.,  v.  63,  pp.  75-80.,  January 
1978.  ~  ~  . . . . 


165 


Commander 

Naval  Security  Group  Command 
Naval  Security  Group  Command  Headquarters 
3801  Nebraska  Avenue,  NW 
(Attn:  Code  G82) 

Washington,  DC  20390 

Professor  W.  Ray  Vincent,  Code  62Ja 
Department  o-P  Electrical  and  Computer 
Eng i neer i ng 

Naval  Postgraduate  School 
Monterey,  CA  93943 

Director,  Research  Administration 
Code  012 

Naval  Postgraduate  School 
Monterey,  California  93943-5000 


